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Abstract 

HiggsBounds is a computer code that tests theoretical predictions of models with ar- 
bitrary Higgs sectors against the exclusion bounds obtained from the Higgs searches 
at LEP and the Tevatron. The included experimental information comprises exclu- 
sion bounds at 95% C.L. on topological cross sections. In order to determine which 
search topology has the highest exclusion power, the program also includes, for each 
topology, information from the experiments on the expected exclusion bound, which 
would have been observed in case of a pure background distribution. Using the pre- 
dictions of the desired model provided by the user as input, HiggsBounds determines 
the most sensitive channel and tests whether the considered parameter point is ex- 
cluded at the 95% C.L. HiggsBounds is available as a Fortran 77 and Fortran 90 
code. The code can be invoked as a command line version, a subroutine version and 
an online version. Examples of exclusion bounds obtained with HiggsBounds are 
discussed for the Standard Model, for a model with a fourth generation of quarks 
and leptons and for the Minimal Supersymmetric Standard Model with and without 
CP-violation. The experimental information on the exclusion bounds currently im- 
plemented in HiggsBounds will be updated as new results from the Higgs searches 
become available. 
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Program Title: HiggsBounds 

Journal Reference: 
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Licensing provisions: 

Program requirements: HiggsBounds can be built with any compatible Fortran 77 
or Fortran 90 compiler. The program has been tested on x86 CPUs running under 
Linux (Ubuntu 8.04) and with the following compilers: The Portland Group Inc. 
Fortran compilers (pgf77, pgf90), the GNU project Fortran compilers (g77, gfor- 
tran). 

Programming language: Fortran 77, Fortran 90 (two code versions are offered) 
RAM: minimum of about 6000 kbytes (dependent on the code version) 
Keywords: Elementary Particle Physics; General, High Energy Physics and Com- 
puting; Higgs bosons; Higgs search; LEP; Tevatron; Beyond the Standard Model. 
PACS: 14.80.Bn; 14.80. Cp; 12.60.Fr. 
Classification: 

External routines/libraries: HiggsBounds requires no external routines/libraries. 
Some sample programs in the distribution require the programs FeynHiggs 2.6.x 
or CPsuperH2 to be installed. 

No. of lines in distributed program, including test data, etc.: 61645 
No. of bytes in distributed program, including test data, etc.: 13721912 
Distribution format: tar.gz 

Nature of problem: Determine whether a parameter point of a given model is ex- 
cluded or allowed by LEP and Tevatron Higgs-boson search results. 
Solution method: The most sensitive channel from LEP and Tevatron searches is de- 
termined and subsequently applied to test this parameter point. The test requires 
as input model predictions for the Higgs-boson masses, branching ratios and ratios 
of production cross sections with respect to reference values. 

Restrictions: In the current version, results from decay-mode independent Higgs 
searches and results of searches for charged Higgs bosons are not taken into ac- 
count. 

Running time: about 0.01 seconds (or less) for one parameter point using one pro- 
cessor of an Intel Core 2 Quad Q6600 CPU at 2.40GHz for sample model scenarios 
with three Higgs bosons. It depends on the complexity of the Higgs sector (e.g. the 
number of Higgs bosons and the number of open decay channels) and on the code 
version. 
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LONG WRITE-UP 



1 Introduction 

A major goal of the particle physics programme at the high energy frontier, 
currently being pursued at the Fermilab Tevatron collider and soon to be 
taken up by the CERN Large Hadron Collider (LHC), is to unravel the na- 
ture of electroweak symmetry breaking. While the existence of the massive 
electroweak gauge bosons (W ± ,Z), together with the successful description 
of their behaviour by non-abelian gauge theory, requires some form of elec- 
troweak symmetry breaking to be present in nature, the underlying dynamics 
is not known yet. An appealing theoretical suggestion for such dynamics is the 
Higgs mechanism [1] , which implies the existence of one or more Higgs bosons 
(depending on the specific model considered). Therefore, the search for Higgs 
bosons is a major cornerstone in the physics programmes of past, present and 
future high energy colliders. 

Many theoretical models employing the Higgs mechanism in order to account 
for electroweak symmetry breaking have been studied in the literature, of 
which the most popular ones are the Standard Model (SM) [2] and the Mini- 
mal Supersymmetric Standard Model (MSSM) [3J. Within the SM, the Higgs 
boson is the last undiscovered particle, whereas the MSSM has a richer Higgs 
sector, containing three neutral and two charged Higgs bosons. Among alter- 
native theoretical models beyond the SM which lead to the prediction of Higgs 
bosons, the most prominent are the Two Higgs Doublet Model (THDM) [I], 
non-minimal supersymmetric extensions of the SM (e.g. extensions of the 
MSSM by an extra singlet superfield [S]), little Higgs models [6] and mod- 
els with more than three spatial dimensions [TJ. Furthermore, several models 
have been proposed which mainly extend the Higgs sector of the SM, e.g. ex- 
tensions of the SM by electroweak singlet scalar fields j8] , or the Private Higgs 
Model [9] which introduces one Higgs doublet for each fermion generation, or 
models with Higgs- dependent Yukawa couplings [TU] . 

LEP has searched for the SM Higgs boson [TT], the MSSM Higgs bosons |12j . 
and other, more "exotic", manifestations of the Higgs mechanism [13]. After 
the termination of LEP, the search has been taken up again by the Tevatron. 
So far, no signals of Higgs bosons have been found, and LEP and Tevatron 
turned the non-observation of Higgs signals into cross section constraints. The 
constraints are provided by experiments in the form of limits on cross sections 
of individual signal topologies (such as e + e~ — > hiZ — > bbZ or pp — > hiZ — > 
bbl + l~) or in the form of combined limits for a specific model, such as the 
SM. In the latter case, the individual topological cross sections have been 
combined using the proportions of the individual contributions as predicted 
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by the model. 



An important part of the phenomenological exploration of physics models be- 
yond the SM is the confrontation of such models with existing experimental 
constraints. In particular, checking the predictions of the model's Higgs sector 
with bounds from experimental searches for a Higgs boson is vital. In order 
to test whether the parameter space of a certain model is compatible with the 
limits from the Higgs searches at LEP and the Tevatron, it is, in general, not 
possible to employ the experimental results that have been obtained for the 
Higgs search in the SM (or any other specific model). This is because, in a gen- 
eral model of new physics, individual Higgs signal topologies will contribute in 
different proportions than, for instance, in the SM. Thus, for models for which 
no explicit experimental analysis has been carried out, or, like in the case of the 
MSSM, for model parameters that differ from the benchmark values chosen in 
the experimental analyses, one needs to resort to the limits on the individual 
topological cross sections provided by the experimental collaborations. These 
limits exist for various search channels from LEP and the Tevatron (the latter 
ones are frequently updated as new data becomes available). Comparing the 
predictions of a particular model with the existing experimental bounds on 
the various search topologies can be quite a tedious task as it involves the 
implementation of experimental results that are distributed over many differ- 
ent publications and combining these results requires a procedure to ensure 
the correct statistical interpretation of the exclusion bounds obtained on the 
parameter space of the model. 

We present here the program HiggsBounds, which is a tool designed to facili- 
tate the above task so that wide classes of models can easily be checked against 
the state-of-the-art results from Higgs searches. This should be useful for ap- 
plications in Higgs phenomenology and model building. HiggsBounds takes 
theoretical Higgs sector predictions, e.g. for a particular parameter scenario of 
a model beyond the SM, as input and determines which Higgs search analysis 
has the highest exclusion power according to a list of expected exclusion limits 
from LEP and the Tevatron (an expected exclusion limit corresponds to the 
bound that one would obtain in the hypothetical case of an observed distribu- 
tion that agrees precisely with the background expectation) . In order to ensure 
the correct statistical interpretation of the obtained exclusion bound as a 95% 
C.L., the comparison of the model with the experimental limits has to be re- 
stricted to the single channel that possesses the highest statistical sensitivity 
for setting an exclusion limit. For this channel, the program then compares 
the theoretical prediction for the Higgs production cross section times decay 
branching ratio with the actual experimental limit and determines whether 
or not the considered parameter point of the model is excluded at 95% C.L. 
So far, no bounds on charged Higgs bosons are taken into account. This will 
be included in a later version. HiggsBounds is model-independent, in that it 
can be used for models with arbitrary Higgs sectors. The input needed is the 
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number of neutral Higgs bosons in the model and masses, decay branching 
ratios and ratios of production cross sections versus corresponding reference 
(usually SM) values for all neutral Higgs bosons in the model. The code has 
both a Fortran 77 and Fortran 90 version. It can be operated in a command 
line mode that can process input files in a variety of formats, as a subroutine 
suitable for inclusion in user applications, and as an online version, available 
at the URL: |http : //www . ippp ■ dur . ac . u k/HiggsBounds , HiggsBounds in- 
cludes sample programs which demonstrate its usage. Sample programs which 
demonstrate how HiggsBounds can be used in conjunction with the widely 
used programs FeynHiggs [Tljl5lll6|[TT] and CPsuperH [TS] for Higgs-sector 
predictions in the MSSM are provided by default. 

The rest of the paper is organised as follows. In Section [2] we describe our 
general approach for the implementation of experimental Higgs constraints 
followed by specific details on the implementation of LEP and Tevatron results. 
In Section [31 the general operating instructions for HiggsBounds are specified 
and accompanied by small sample programs, while in Section H] some examples 
of use for HiggsBounds are described. Section contains a summary and 
outlook and is followed by an appendix which contains details on the definition 
and statistical interpretation of the confidence limits. 



2 Implementation of experimental Higgs constraints 

2.1 General approach 

Tevatron and LEP turn(ed) the non-observation of Higgs signals into 95% C.L. 
upper limits on cross sections for individual signal topologies. These limits are 
given as functions of the anticipated mass of the Higgs boson. By a certain 
Higgs signal topology X, we mean a specific combination of a single (double) 
Higgs boson production process, P{h) (P(hi, ^2)), and decay final state(s) F 
of the Higgs boson(s) h (hi, h-2)- In the limit of a narrow-width Higgs boson the 
cross section a(X) of the signal topology X factorises into the on-shell Higgs 
production cross sections u(P) times the appropriate Higgs decay branching 
ratio(s) BR(h -> F) (BR(h u h 2 -> F)). 

The Tevatron and LEP cross section limits are understood to be applicable to 
models which do not change the signature of the background processes con- 
siderabljQ which do not significantly change the kinematical distributions of 



This is not such a strong extra restriction on new physics models as it may seem, 
because viable models are required to agree with present experimental results which 
are well described by the SM, such as the LEP and SLD EW precision data [19] and 
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the signal topology X (e.g. rj, px distributions of the final state particles) from 
what has been assumed in the corresponding analysis. This is an unavoidable 
limitation of the standard approach adopted in the literature of expressing 
search results in terms of limits on total cross sections times branching ra- 
tios for certain signal topologies. Detection efficiencies, for instance the b-jet 
tagging efficiency, usually depend on the kinematical configuration of the sig- 
nal event (particularly the r\ and px of the final state particles). Therefore, a 
dedicated experimental analysis for a model giving rise to kinematical distri- 
butions that strongly differ from the ones assumed in the analysis on which 
the published cross section limit is based would lead to a somewhat modified 
cross section limit. For instance, applying the cross section limit of a Tevatron 
analysis which relies on b-jet tagging to a model where the distribution of 
the signal b-jets differs strongly from the one assumed in the Tevatron anal- 
ysis would lead to a slightly too weak/strong limit if a higher/lower fraction 
of the signal b-jets than assumed in the analysis is radiated into the central 
part of the detector where the b-jet tagging efficiency is higher than in re- 
gions closer in angle to the beam axis. HiggsBounds provides three options 
for passing model predictions to the program: (hadr, part and effC). The 
extent to which the kinematical distributions may vary from those assumed 
in the analyses depends on which of these input options is chosen. We will 
discuss this issue in further detail in the sections describing the input options. 

The Tevatron and LEP cross section limits are usually given for a narrow- 
width Higgs boson. From the narrow-width limits, some experimental analyses 
extract limits for Higgs bosons with a substantial decay width by representing 
the signal of a "large-width" Higgs boson as a weighted sum of mock narrow- 
width Higgs signals. Examples of this procedure can be found in [2~ni22ll23|T2"l] . 
In particular, [23] shows that the effect of a non-vanishing width on the ex- 
tracted limit can be parametrised by a function of the Higgs mass rrih and 
the ratio T tot {h)/mh- If model-independent information of this kind also be- 
came available (for a fine grid of points) for other analyses, this would greatly 
facilitate the task to implement width-dependent Higgs exclusion limits. 

At present, we constrain the implementation of HiggsBounds to the case of 
narrow- width Higgs boson exclusion limits. Neglecting the width of a Higgs bo- 
son usually leads to tighter bounds on the topological cross sections than those 
obtained through taking the width into account. The user should be aware of 
this fact when considering exclusion bounds on parameter regions that give 
rise to very large Higgs-boson widths. In a future version of HiggsBounds we 



the precision measurements of gauge boson production at LEP and the Tevatron 
[20] . These analyses are much more sensitive to changes in the background to Higgs 
searches than the Higgs searches themselves. Therefore, the models with new physics 
in the Higgs sector which have been considered in the literature usually do not show 
strong deviations from the SM in the background processes. 
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plan to incorporate a proper width treatment for cases where width-dependent 
results are provided by the experimental collaborations. As mentioned above, 
HiggsBounds is not meant to be applied to models for which the background 
rates to the considered signal processes differ drastically to the SM case. 

For some specific models (the SM and some benchmark scenarios of the 
MSSM) the results of several Higgs signal topologies have been combined by 
the experimental collaborations according to the relative contributions to the 
total cross section. In the combination, a detailed knowledge of the overlap 
between the individual experimental searches is used, which is not available 
using the limits on topological cross sections only. Hence, a combined limit, for 
instance for the SM, is more sensitive than a combination of individual topo- 
logical limits. The use of detailed SM dependent information during the cal- 
culation of the combined limit prohibits its application to those Higgs bosons 
of new physics models which are not predicted to have all relevant couplings 
proportional to the SM values. However, such models can still contain one 
(or some) Higgs boson(s) which are SM-like. For those Higgs bosons, the SM 
combined limit can be applied. 

The basic flow of HiggsBounds is as follows. 

(1) The user provides (in a form to be specified below) for all neutral Higgs 
bosons hi (i — 1, . . . ,nHi ggs ) of the considered model the predictions for: 
Higgs masses, total decay widths, the relevant branching ratios for the 
decay of Higgs bosons hi into known ordinary particles (OP, i.e. SM 
fermions and gauge bosons) and into pairs of other Higgs bosons hj, and 
the relevant cross sections for Higgs production processes P normalised 
to reference (usually SM) values: 



The precise definitions of the reference cross sections are given in the LEP- 
and Tevatron-specific subsections below. The user can choose to specify 
all the above quantities directly, or, in the simplest form, one just needs 
to provide m^, T tot (hi), BR(/ij — > hjhj) and effective couplings of the 
Higgs bosons to ordinary particles. The remaining branching ratios and 
cross section ratios listed above can then be calculated from this input in 
the effective coupling approximation. For the preparation of normalised 
input, the user may apply the SM predictions for Higgs boson production 
cross sections and decay branching ratios which HiggsBounds provides 
as external routines (see Table [T2l in Section [373]) . 
(2) With the above model input, the program calculates the model predic- 
tions for the cross section (normalised to a reference value in some cases) 





a rci (P(H),m H ) mH=mh ' a Te {(P(H, H'),mH:rn H >) 
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for each Higgs signal topology Q m odei(X). Depending on whether the 
exclusion result for a particular search topology has been given by the 
experimental collaborations as a relative or absolute limit, the program 
evaluates 

[ax BR] mode i r , 

V model — — f fpp-j or F x DtX mo del > 

[a x BR] ref 

where a and BR denote the production process and decay branching ratio 
of the Higgs boson, respectively. 

(3) From the experimental results, we read off the values Q bs(X) and Qex P ec(AT) 
corresponding to the observed and expected 95% C.L. exclusion bound on 
the quantity Q(X), respectively. To that end, we have implemented many 
of the available observed and expected exclusion limits (preliminary and 
final results) from LEP and the Tevatron as data tables which are read 
in by the program during start-up. In order to provide values for Q b s (X) 
and Qexpcc(^) for continuous Higgs mass values, the program interpolates 
Q-values linearly between neighbouring Higgs mass points given by the 
experimental collaborations. Details on the definition of 95% C.L. limits 
can be found in Appendix [A] 

(4) The program determines the search topology Xq with the highest ratio of 
Qmodei to expected Qexpec, i-e. out of all topologies X, it finds the channel 
X where 

Qmodel(^) 
Qexpec (-^0 

is maximal. This topology has the highest statistical sensitivity for ex- 
cluding the model prediction. 

(5) We then look at the ratio of Q mo dci to the observed Q Q b s value for this 
channel X and, if 

Qmodcl(^o) > j ^ 
Qobs(^o) 

this particular parameter point is excluded by the corresponding experi- 
mental analysis for the search channel Xq at 95 % C.L. In order to ensure 
the correct statistical interpretation of the obtained exclusion bound as a 
95% C.L., it is crucial to perform the test of Eq. (j2J) only for exactly one 
search channel (namely the one with the highest statistical sensitivity), 
and not for all available search channels. In the latter case the derived 
constraint would in general not correspond to a constraint at 95 % C.L., 
due to the maximally 5 % probability of each individual comparison of 
Qmodei and Qobs to yield a false exclusion. 

As additional information, the web version of the program HiggsBounds also 
informs the user which topologies had the second and third highest statistical 
sensitivity and whether or not the ratio of Q mo dei to observed 95% C.L. limit 
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Qobs w as greater than one for those topologies. However, this information does 
not contribute to the decision about whether or not this parameter point is 
excluded. 



There are a few more general implementation details which should be men- 
tioned here. 

Using SM combined limits, SM-likeness 

In order to use the best available exclusion limit, it makes sense to use the 
results of the SM combined analyses for any Higgs boson predicted to behave 
in a similar way to the SM Higgs boson (in terms of its production and decay). 
A SM analysis, combining several signal topologies assumes that the relative 
contribution to the event rate by each included search topology is according to 
the proportions in the SM. For a given Higgs mass m#, the only free parameter 
left in such analyses is one overall scale factor, which multiplies all Higgs signal 
cross sections and which is varied in order to determine the value for which 
95% CL. exclusion occurs. 

Therefore, the results of a SM analysis, combining M signal topologies, are 
exactly applicable to a Higgs boson h of an alternative model if for all of the 
M search topologies X n the ratio of topological cross sections (model versus 
SM) is independent of X n , i.e. 

"model ( n _ cong |-_ = ■ a ^ n = 1, . . . , M , (3) 

cr SM(^n) 

with a proportionality constant a, or equivalently 

"model "SM (Xn) 



""model 

(X k ) o- SM (X k ) 



for all n, k e {!,..., M} . (4) 



For a model which passes the conditions (jSJ) or (jlj) to sufficient accuracy for a 
set of M Higgs search topologies, the cross section prediction for a combination 
of those topologies, normalised to the SM reference value, fulfils 



Xm=l "model (X r , 
J2n=l CsM(X n ) 



r\ (iv V ~\ \ — ^n=l "model \^n) _ 

VmodellVM) • • • A M/J •— ^ M ~ — Ct . [0) 



Here and in the following, it is always understood that the mass of the SM 
Higgs boson m# is set to nif L . In the narrow-width approximation, which 
we assume to hold throughout (see the discussion above), each cross section 
cr(X n ) of the Higgs signal topologies X n factorises into a production cross 
section a(Pi(h)) and decay branching ratio BR(/i — > F k ), where % and k are 
functions of n, and n — 1, . . . , M. Note that the index n labels distinct signal 
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topologies, i.e. specific combinations of one Higgs production process with one 
decay final state. SM-likeness is fulfilled if 



a(Pi(h)) 

model BR mo dei(^ - ► F k ) 
oiPjHjU BR SM (H^F k ) ~ C ° nSt - W 

holds to sufficient accuracy for all M topologies. If this is the case for a partic- 
ular Higgs boson of a certain model, a SM analysis combining precisely those 
M channels also applies to that Higgs boson. 

In order to decide whether or not a Higgs boson of a certain model is suf- 
ficiently SM-like such that the corresponding SM combined analysis can be 
applied, we use in HiggsBounds a pragmatic definition of SM-likeness. For 
each of the N cs distinct production cross sections cT modcl (Pj (/?,)) and the N BR 
distinct decay branching ratios BR(/i — > F k ) which appear in the list of M 
signal topologies, we determine the normalised mean value and the deviation 
from the mean: 



I Ncs 









b = 






5b k = 


h-b, 




bk = 


BR mo dci(^ - 


+ F k ) 


BRsm(-^ — 


>F k ) 



Ssi Si s , 

Wlths ' = ^.W(H)) ' i= 55 • ( ) 

The approximate SM-likeness of the given model is then ensured if the devi- 
ations from a constant expression according to Eq. (jBJ) are small. For a given 
combination (i, k) of a production and decay process, the left hand side of 
Eq. (|5J) can be expressed as: 

Qmodci(PiQ)) BR modc i(fe -> F k ) _ - ( 5s,j 5b k 5sj5b k \ 

a SM (P t (H)) BR SM (H^F k ) s ^ s0 {^ g + ~ b + 5 g J ■ 

(10) 

In the program, we consider a corresponding SM analysis applicable if the 
maximum relative deviation from the mean value sb satisfies 

/ Ssi 5b k 5si5bk\ , s 

max - + ^ + ^ <e, 11 
i,k \ s b sb J 

where we have chosen e = 2%, i.e. the predictions for the different topolog- 
ical cross sections, normalised to the SM values, are required to be equal to 
each other with at least 2% accuracy. This restriction is quite stringent and 
may cause a SM combined analysis to be not applied to some models where 
it may actually be justifiable. The user can remedy this situation by chang- 
ing the setting of the corresponding variable eps in the internal subroutine 
check_SM_likeness (Fortran 77) or SMlikeness (Fortran 90) in the code. 
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So far, HiggsBounds uses analyses with SM combinations of topologies only 
from the Tevatron experiments. This is the most important CcLSGj clS CDF and 
D0 have not done combinations of their results for individual search topologies 
so far, while the LEP results we are using represent the combination of results 
from all four LEP experiments for each individual channel. 

A more elaborate criterion for determining whether SM combined analyses 
are applicable to Higgs bosons of models beyond the SM can be sketched as 
follows. The variation of the ratios of signal topology cross sections (model 
versus SM) due to non-exact SM-likeness of the model prediction should stay 
within a fraction of the la band of the sensitivity of the observed limit of the 
analysis with respect to signal-like fluctuations. The fixed maximum deviation 
of 2%, described above, lies well within the typical expected statistical fluctu- 
ation of the limits in the SM combined analyses. We plan to implement such 
a more elaborate criterion in future updates of the program. 

Adding of signal channels 

In models with more than one Higgs boson it can happen that some of the 
Higgs bosons, hi, have masses which are quite close to each other, say 

\m hi - m hj \ < Am h (12) 

for certain i and j. This is the case for instance in the parameter space of the 
(CP-conserving) MSSM for the heavy Higgs bosons H° and A for Ma 3> Mz- 
If the masses are so close to each other that their signals appear within a cer- 
tain mass window, it may happen that the two signals are not experimentally 
distinguishable anymore, or, at least, that it makes sense to combine their 
signal contribution into one signal to be compared with the data. 

In the narrow-width approximation (which is used in HiggsBounds, see the 
discussion above) it is assumed that the total widths of the Higgs bosons 
Ttot(^i) are all (much) smaller than their mutual mass differences, 

r to t(^i) < \m hi - m hk \ for all % and k^i . (13) 

In this case the signal channels factorise into production cross section and 
decay branching ratio, so that the addition of different signal cross sections 
may formally always be carried out where it is necessary. In parameter regions 
where the narrow-width approximation does not hold, however, one may worry 
about the possible effects of interference terms arising from squaring the am- 
plitude of the full process involving both production and decay. It should be 
noted that even in cases where the total widths of two Higgs bosons are sig- 
nificantly larger than their mass difference it can happen that interference 
terms are small. An example for such a situation is the above-mentioned case 
of the MSSM with real soft-breaking parameters for Ma S> Mz- While for 
large tan (3 the widths of the two Higgs bosons can be very large, the different 
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quantum numbers of the CP-even and the CP-odd state nevertheless forbid 
an interference contribution. 

In determining the cross sections o~(X) for a given signal topology X, Higgs- 
Bounds does the following. Suppose there are n mggs neutral Higgs bosons hi 
in a given model and let the numbering be according to a mass order, i.e. 

m hl < m h2 < ■ ■ ■ < m hnmggs . (14) 

Then, HiggsBounds works through a loop from i — 1 to ^mggs and adds up 
al0 production cross sections of Higgs bosons hj (j > i) which satisfy 

m hj - m hi < Am h (X) . (15) 

In the implementation, the X-dependence of Am^ is limited to two different 
values: one for LEP and one for Tevatron search topologies, represented by 
the variables delta_Mh_LEP and delta_Mh_TEV, respectively. Adding of Higgs 
signal cross sections is mainly relevant for Tevatron Higgs search analyses, 
where the invariant mass distributions of the potential Higgs boson decay 
products usually use bin sizes of the order of 10 GeV. For the LEP analyses, 
the Higgs mass resolution is typically around 2-3 GeV |llj . For the LEP results 
the addition of Higgs signal cross sections is implemented in HiggsBounds 
only for the Higgsstrahlung process where the Higgs boson decays directly 
into SM particles. As default values for delta_Mh_LEP and delta_Mh_TEV we 
have chosen (in units of GeV): 

delta_Mh_TEV = 10 , 
delta_Mh_LEP = 2 . 

The default setting of delta_Mh_TEV for the Tevatron search channels has also 
been tested against known MSSM Higgs search results [23|24II25] which we 
could reproduce successfully with HiggsBounds within the limitations of the 
narrow-width approach. 

In the calculation of normalised quantities, the reference cross section and 
branching ratio are evaluated at the average mass of the Higgs bosons con- 
tributing to a signal topology. 

Provided SM normalisation 

In the evaluation of the model predictions Q mo de\(X) and <5modci({Vi, . . . , Xm}) 
for single and multiple search topologies, respectively, to be compared with 
the expected and observed limits Q b s and Qexpec, HiggsBounds uses inter- 
nally SM predictions for Higgs production cross sections and decay branching 

2 If the cross section for a SM combination of signal topologies is calculated, pro- 
duction cross sections of those Higgs bosons are only included in the sum if they 
also fulfil the SM-likeness criterion. 
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ratios. Out of the set of implemented Higgs search results, the need for us- 
ing such predictions internally arises only for Tevatron analyses. The provided 
SM predictions for Higgs boson production cross sections and decay branching 
ratios use results of HDECAY 3.303 [25] and the TEV4LHC Higgs Working 
Group [27] (see Table [T2l below for references to the original publications). In 
addition, for the bH production cross section we need several predictions with 
different kinematic cuts on the b jet which are not available from [27]. We 
used HJET 1.1 [28J to calculate those, but checked that the uncut cross sec- 
tion closely resembles the result reported in [27] (see Section I3~31 for details). 
If the model input is given in terms of ratios of partonic cross sections, the 
calculation of hadronic Higgs production cross sections is facilitated by using 
internally SM cross section ratios which we have calculated (see Section 12.3.21 
for details). If the model input is given in terms of effective couplings, we make 
use of results from VBFNLO [2H] in order to facilitate the calculation of the 
Higgs production cross section via vector boson fusion (see Section 12.3.31 for 
details). 

The rationale behind the choice of SM normalisation is that virtually all Teva- 
tron analyses implemented in this program use these predictions when normal- 
ising their cross section limits to SM quantities. Thus, describing deviations 
from the SM of a new model by using the SM normalisation the experimen- 
tal analyses have chosen allows for the most accurate interpretation of the 
limits. By using SM predictions for Higgs production cross sections which de- 
viate from the internally used ones by a certain percentage (be it because of 
a different loop order, different numerical values of input parameters, renor- 
malisation scheme or choice of parton distribution functions), a deviation of 
the same relative size will be caused in the quantities Q model- The user should 
bear this in mind when interpreting the output of HiggsBounds. 

2.2 LEP limits 

The main Higgs boson production processes at LEP have been Higgsstrahlung 
and double Higgs production: 

e + e~ — > hkZ , e + e~ — > hkK . (16) 

Both processes are assumed to be mediated via s-channel Z boson exchange 
in leading order of perturbation theory. 

Currently, HiggsBounds uses the exclusion limits on topological cross sections 
obtained from the combination of all LEP2 and some LEP1 Higgs boson search 
results from Ref. [12] which considers the Higgs decay final states bb and 
t + t - . Apart from SM search topologies, this analysis also considers double 
Higgs production topologies, and almost all possibilities for the decay of an 



13 



search topology X 


reference 






— > \00)Zj 




11 OIQfll 

[izpoUj 


f r,k1 Q 1 /IK 




^ ('A \7 


-> (T T T 




11 OIQfil 

|lz|oU] 


table 14c 


e~*~e 


_v / A, 

->• (Aifc — > 


hihi)Zj — 


[0000 )Zj 


11 OIQ1 1 
|±Z|ol] 


table lo 


e~*~e 


_v / A, 

->• [h k — > 


h h \ 7 

hihi)Zj — 




11 OIQ1 1 

[izjol] 


table id 


e~*~e 


-> {hktii) 


— > (J3000) 




11 OIQ1 1 

|lz|ol] 


table lo 




-> {tlklli) 


— ► (T T 


T T J 


11 OIQ1 1 

|lz|ol| 


table iy 


6 6 




hh\h- 
III III Jill — 


r \JJOOO JOO 


1 9 Q 1 


fnUlo on 

table zu 


e + e~ - 


-> (/i fe -> 


hihi)hi - 


> (t + t~t + t~)t + t~ 


|12|31| 


table 21 


e + e~ 


->• (/ifc -> 


hihi)Z — 


(66)(r+r-)Z 


mi 




e + e~ 




bb)(h t - 


T+T") 


m 




e + e~ 


-> (/i fc -> 


r+r-)(fti 


66) 


m 





Table 1 

LEP search topologies used by HiggsBounds. For each Higgs search topology which 
has been studied in Ji2| . we indicate the corresponding table number. It is assumed 
in all tables that m/ lk > . . In the program, we use more fine-grained results \3C\31\ 
from the analysis Ji2f than have been displayed as tables in \12jj . We also use equally 
fine-grained results for the expected limits \30\31\ . which have only been displayed 
in graphs in \12\j . and for the observed and expected limits of a few channels, which 
were not published \31}j . 

on-shell Higgs boson into two lighter one and their subsequent decay into 
matter particles. Because of the limitation in [12] to bb and r + r~ final states, 
Higgs bosons which decay predominantly invisibly cannot be well constrained 
by HiggsBounds at present. We will include the appropriate LEP results in 
future updates of the code. The search topologies included in HiggsBounds 
up to now are shown in Table [TJ 

The limits on topological cross sections in [12] have been given in the nor- 
malisation where all Higgs decay branching ratios are equal to 1. For given 
Higgs mass(es), limits on the production cross sections of the individual search 
topologies X have been derived for this normalisation in the form of scaling 



3 The only possibilities not analysed in [12] are the ones where more than two 
distinct Higgs bosons appear. For example, this occurs if a primarily produced Higgs 
boson hk decays into two different Higgs bosons hj and hi, or if two Higgs bosons 
hkhj are primarily produced and one of these, say hk, decays into two Higgs bosons 
hi k all different). We are not aware of any studies by the LEP experiments of 
those more general topologies. 
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factors Sq5(X), defined as 



0"n 



.XX)/a ief (X), 



(17) 



where <J max is the largest cross section compatible with the data at 95 % C.L. 
and <7 re f is a reference cross section for the Higgs production process. The 
quantities S 95 (X) of [T2] correspond to the quantities Q(X) mentioned in the 
description of the program flow in Section 12. 1[ 

For Higgsstrahlung processes, e + e~ — > h^Z — > (final state)Z, the reference 
cross section <7 re f is the SM cross section 



cr rcf (Higgsstrahlung) = a(e + e~ -> HZ) SM (18) 
with the SM Higgs mass m# chosen as mh k - 

For double Higgs production processes, e + e~ — > H'H — > (final state), there 
is no direct SM reference process. However, as in [12], we choose a reference 
cross section for a fictitious production process of two scalar particles (H ; , 
H) with masses mw = fnh k and m# = via a virtual Z exchange with a 
standardised coupling constant 



9lPHZ 



(19) 



where e denotes the electromagnetic coupling constant, and s w and cv the sine 
and cosine of the electroweak mixing angle, respectively. The cross section of 
this process in leading order is then completely determined by the Higgs masses 
and SM input and related to the SM Higgsstrahlung cross section via a simple 
phase space factor: 



a Te f(H' H production) = A (m H /,m H , s) &Hz( m H) 



\(m H >,m H ,s) 



X]i 2 z (s [s) + l 

2i 



(20) 



1 - 



1 - 



2i 



This reference cross section coincides with the MSSM or THDM tree-level 
cross section for the process e + e~ — > h°A° if the Higgs mixing-angle dependent 
factor cos(/3 — a) is divided out and m^o and m^o are chosen as and m#. 

In order to apply the LEP Higgs constraints implemented in HiggsBounds 
to a model, the user has to provide model information which allows model 
predictions to be calculated for all Higgs search topologies X listed in Table 
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m 



r\ / -\/-\ ^modcl 

(P(X))BR model (F(X)) 

WmodellAJ - , [ZL) 

a iei {P{X)) 



where the symbol P(X) denotes the corresponding Higgs production process 
and F(X) the corresponding Higgs decay chain, e.g. h k — > bb or h k — > hihi —> 
bbr + r~ . The quantity Q mo dei (X) is compared with the experimental limits 
Qexpec(^) and Qobs(-^) according to the description in Section |2~T1 Currently 
HiggsBounds allows model information to be provided in three ways. 



2.2.1 LEP input (HiggsBounds Input Option part and hadr^): cross section 
and branching ratios 

In order to be able to apply the implemented LEP limits to the model under 
study, the user is asked to provide model predictions for the Higgs boson 
masses (in units of GeV), for the branching ratios 

BR mo dei(^i — ► bb) , BR modc i(/ij — »■ t + t~) , BR mode i(/i fc — > hihi) , (22) 

and for the ratios of cross sections 



o"modei(e + e — > h k Z) o"modei(e + e — > h k hi) 



a ie{ (e+e~ HZ) ' cr re f(e+e _ — ► H'H) 

for k,i G {1, . . • ,n ffig gs}- 



(23) 



With this information, HiggsBounds obtains model predictions Q m odei(^0 f° r 
all topologies X specified in Table [TJ For example, the Q mo dci value for the 
channel e + e~ — > h\Z — > (bb)Z is given by (suppressing the initial state in the 
notation): 

Qmodci Uh x )Z -> (bb)z) = ffmodc '^f BR modcl (hi - 66) , 
v ' a rc{ {HZ) 

with ran — m h 1 i an d similarly for the channel e + e~ — > (h 2 — > h\h\)Z — > 
(bbbb)Z it is: 

Qmodci ((h 2 -> /ii/ii)^ -> (bbbb)Z 



^BR modcl (/i 2 — > /ii/ii) (BR modcl (/i! — > bb)^ 



Vref(HZ) 

with m# = i7ih 2 . 

The HiggsBounds Input Option hadr amounts to exactly the same LEP re- 
lated input as described above. It differs for the Tevatron related input (see 
Section [2.31) . 
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If the model under study predicts kinematical distributions which differ sig- 
nificantly from the distributions of the reference processes (1180 or (|20|) . the 
exclusion result, returned by HiggsBounds, can only be considered an esti- 
mate0 



2.2.2 LEP input (HiggsBounds Input Option effC): effective couplings 



The user is asked to provide model predictions for: 



m hk , T tot (h k ) 



(„model\ 2 /„modcl\ 2 / model 

9h k zz \ I 9h k hjZ \ I ^hj, 

dWzz ) \9whzJ \ 9 H j 



h k ff,eS 



,SM 
J Hff 

BR mo dci(^fe — ► KK) , (24) 



for k,i G {1, . . . , 77-Higgs} and / G {b, r}. The normalised effective squared 
coupling of the Higgs boson fermion interaction is defined as an approximation 
for the ratio of partial decay widths, 
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model 
h k ff,cS 

SM_ 
y Hff 



"pmodel I ^y-, 



H^ff 



(m H ) 



(25) 



m H =m hk 

and the reference coupling constants are defined by Eq. ( fT9l) and by 

2 



(9 



HZZ) 



(9 



HffJ 



e mw 

1 e rrif V 

2 s w m w J 



(26) 
(27) 



where my/ and rrif denote the masses of the W boson and fermion /, respec- 
tively 



For the purpose of running HiggsBounds, only the ratio of partial widths cor- 
responding to the right-hand-side of Eq. (1251) is needed. In order to completely 
specify effective couplings for neutral Higgs boson interactions with fermions, 
in general a scalar and a pseudoscalar part of the coupling is needed. The 
Feynman rule for the coupling of a generic neutral Higgs boson h to fermions 
can be written as: 



G(hff) = i{g s l + Z&75) , (28) 



In such a case, the general cross section limits given by the experimental col- 
laborations would have to be replaced by a dedicated experimental analysis of a 
particular model giving rise to these kinematical distributions, taking into account 
detailed information about the detector response to the signal events. 
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where g s and g p are real- valued coupling constants, and 1 and 75 are the usual 
matrices in Dirac space. A scalar particle, like the SM Higgs boson, has g p = 
and a pseudoscalar particle has g s = 0. If the user has the couplings g s and g p 
available, the user may calculate the normalised effective squared couplings, 
specified in Eq. fl25l) . by using the relation: 



»model 



9 h k ff,eS 
SM_ 
y Hff 



9 



model_ 
s,h k ff 



SM_ 
y Hff 



9 



modcl_ 
P,hkff 



SM_ 
y Hff 



1 



ft 



(29) 



with 



Pf(m hl 



\ 



ml. 



The extra factor of /L 2 in Eq. 
particle it is proportional to 



appears, because the partial width for a 
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scalar particle decaying into // is proportional to /??, while for a pseudoscalar 



The user may also avoid the explicit use of the reference coupling, defined in 
Eq. fl26|) . by using the relation 



( model \ ^ 
9h k ZZ \ 



■p model 
h k -*Zl*)Z 



(*)( m h k 



(30) 



m H =m h 



where the asterisk in brackets indicates the possible off-shellness of the de- 
cay products. Knowledge of the off-shell decay matrix elements allows to use 
relation (|3"U|) for masses smaller than 2Mz- 



From the information specified in Eq. pi}, HiggsBounds calculates the LEP 
part of the cross section and branching ratios required in Input Option part 
in the effective coupling approximation: 



0- m odei(e + e —> h k Z) 



(7 ref (e+e- 

O"model(e + e" 



HZ) 
> h k hi 



h k h 



k'H) 



(„model \ ^ 
9h k ZZ \ 
gTzz) 

( model \ ^ 
9H'HZ \ 
9h'hzJ 



BR mo dei(/ifc — > //) = BR SM (i7 



ff)(m H ) T ^ {mH) 



Ttot(^A 



(31) 
(32) 
(33) 



m H =m hk 



X 



„model 
y h k ff,cS 
SM_ 
y Hff 



for / e {b,r}. 
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From this point, the evaluation of the quantities Q m odc\(X) can proceed exactly 
as exemplified in Section f2. 2. 1[ 

As far as the implemented LEP constraints are concerned, effective couplings 
just rescale all kinematical distributions of the reference processes (ITS]) or (1211 
by common factors. Thus, for the class of models where the Higgs sector can be 
faithfully parametrised in this way, the kinematical distributions correspond 
to those of the reference processes. Consequently, limits on topological cross 
sections can exactly be applied. 



2.3 Tevatron limits 

The main SM Higgs boson production processes at the proton-anti-proton 
collider Tevatron at Fermilab are: 

• Gluon fusion: pp — > H + • ■ ■ via gg — > H, 

• Higgsstrahlung: pp — > HV + ■ • • via qq' — > VH (V = Z, W^), 

• Vector boson fusion (VBF): pp — > Hjj + - ■ ■ via qq' — > V*V'*q"q"' — > Hq"q"', 

• Higgs production associated with heavy quarks: pp — > Hti (Hbb) + • • • , 

where the ellipses indicate the hadronic remainder of the reaction which we 
will suppress in our notation below. 

Currently, the Tevatron experiments CDF and D0 provide separate limits on 
cross sections of individual Higgs search topologies X and limits on combina- 
tions of search topologies according to SM proportions. In particular, combina- 
tions of Higgs search results from both experiments are, so far, only available 
for a SM-like Higgs boson. The results implemented in HiggsBounds corre- 
spond to the Higgs search topologies that have been analysed by CDF and 
D0, using Tevatron Run II data, listed in Tables [2] and [31 New results pro- 
vided by CDF and D0 will be implemented once they appear, thus keeping 
HiggsBounds up-to-date. 

In the absence of a Higgs signal, the Tevatron experiments obtain cross section 
limits from their measurements: <J mSuX (X) is the highest cross section for a 
given assumed Higgs mass mn which is compatible with the only-background 
hypothesis at 95% C.L. Therefore, if a particular model of the Higgs sector 
predicts a cross section <J mo dei(V) for a Higgs boson with mass ran, which is 
higher than the limit, i.e. 

Omodci(A) > a max (X) , (34) 

then the probability that such a prediction is consistent with observation is 
less than 5% and the model parameter point will be called excluded with (at 
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Table 2 

Tevatron analyses for single search topologies, the results of which are used by 
HiggsBounds. The leptons I and I' can be an electron or a muon. References marked 
by an asterisk (*) refer to analyses which have been published or submitted for 
publication. 



least) 95% CL. While some Tevatron Higgs search results are presented in 
the form of absolute cross section limits, a max (X), most of the results are 
presented as limits on signal cross sections normalised to the corresponding 
SM cross section, o"sm(^0H] In the latter case, the exclusion condition ( 1341) 
becomes: 

Cmodel(^) ^ 0"max(^) 

o~sm{X) o- S m{X) 



In order to test the exclusion of a given model by CDF or D0 results according 
to the procedure outlined in Section 12.11 the user needs to provide model 



So far, no Tevatron Higgs search results have appeared which would require to 
use a normalisation to non-SM cross sections. 
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search topologies {X} (analysis) 



reference 



m 
m 



pp -» WH/ZH -^bb + E™ iss - (CDF with 2.3 ft" 1 ) 
pp -> WH/ZH ^bb + E™ 1SS - (D0 with 2.1 ft.- 1 ) 

-» H/HW/HZ/H via VBF,# -» r+r~ (CDF with 2.0 ft." 1 ) 
Combined SM analysis (CDF & D0 with 0.9 - 1.9 ftr 1 ) 
Combined SM analysis (CDF & D0 with 1.0 - 2.4 ft" 1 ) 
Combined SM analysis (CDF & D0 with 3.0 ftr 1 ) 
Table 3 

CDF and D0 analyses which are contained in HiggsBounds and which combine 
search topologies according to SM proportions. 

predictions which allow, for each relevant Higgs boson search topology X, 
either absolute cross sections, 

C^modcl 

(P(X))BR model (F(X)) , (36) 

or cross section ratios, 

r\ I v\ ""model 

(P(X))BR model (F(X)) 
Wm ° dcl1 ] ~ a SM (P(X))BR SM (F(X)) ' [6n 

of hadronic processes to be calculated, where the symbol P(X) denotes the 
corresponding Higgs production process and F(X) the corresponding Higgs 
decay chain. The program allows both variants to be calculated from the 
same input, as the SM predictions for Higgs production processes and decay 
branching ratios are internally available. In Eqs. (1361) and (1371) . we have already 
used the implicit assumption that Higgs bosons have narrow width, for which 
the CDF and D0 limits we employ are valid. We should also stress again that 
the limits are based on analyses which assume that the backgrounds to the 
signal process under study are similar to the backgrounds predicted by the 
SM. For that reason, the branching ratios for the decays of the or the Z 
for a particular individual Higgs search topology X are assumed to cancel in 
the ratio ([3T 



For the calculation of signal cross section predictions where several Higgs 
search topologies {Xi,...,X M } are combined (all cases listed in Tabled]), 
the definition in Eq. (jSJ) applies. In the narrow-width approximation, this 
definition reads: 



2~^n=l "model 

(P(X n ))BR model (F(X n )) 

Zn=l°SM(P(X n )BR SM (F(X n )) 



n> f ( V V 1\ £-in=\ u models \ yl n) ) AJA^modeU/ y^nJJ / Qc 
VmodelU A l5 ' ' ' A M/J — ^JJ ~ , ~, ~ J^TTv \\ ' 
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The quantities Q mo dei are to be compared with the experimental limits Q ex pec 
and Q bs according to the description in Section 12.11 There are three options, 
which allow the user to provide HiggsBounds with the necessary model infor- 
mation in order to calculate the quantities Q mo dei- These are described in the 
following. 



2.3.1 Tevatron input (HiggsBounds Input Option hadr): hadronic cross sec- 
tion and branching ratios 

In order to make use of all the implemented Tevatron limits, the user is asked 
to provide model predictions for the Higgs boson masses m/ lfc (in units of GeV), 
for the branching ratios 

BR mo dcl(^fc —> bb) , BR mo del(/ife —> T + T~) , 

BR modcl (h k -> W+W~) , BR modcl (/i fc -> 77) , (39) 

and for the ratios of the hadronic cross sections 

Qmodei(pp -» h k Z) o- mo dci(pp -> hW^ a mo dei(pP -> h k via VBF) 
&sm(pP -> HZ) ' a SM (pp -> HW*) ' a SM (pp -> H via VBF) 

Cmode\(pP -> h k ) CTmodelipP ~» h k b) 

o- SM (pp^H) ' o- SM {pp^Hb) 
for he {I,.. . ,n ffi ggs}- 

With this information and the internally available SM predictions for Higgs 
production cross sections and decay branching ratios, HiggsBounds calculates 
all values Q mo de\, according to Eqs. ( 1361) . ( |37|) or ( |38|) . which are needed in 
order to test the chosen model against the implemented Tevatron results. 
Note that the branching ratios BR modcl (/i fc — > hihi) are currently not used in 
comparisons with Tevatron results, because, so far, no Tevatron analyses have 
appeared which consider such decays. 

While hadr is the most widely applicable input option, it can also be the one 
that is most difficult for the user to provide, as this option requires the user 
to provide predictions for hadronic cross sections, which involve a convolution 
with parton distribution functions. Therefore, the program offers the Input 
Option part, which requires only ratios of partonic cross sections. This should 
in general be more convenient for the user. However, this input option is 
applicable to a narrower class of models than Input Option hadr. 

If the model under study predicts kinematical distributions which differ sig- 
nificantly from the distributions of the reference SM processes, the exclusion 
result, returned by HiggsBounds, can only be considered an estimate (as dis- 
cussed at the end of Section I2.2.ip . 
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2.3.2 Tevatron input (HiggsBounds Input Option part,): partonic cross sec- 
tion and branching ratios 

The user is asked to provide model predictions for the Higgs boson masses 
rrih k (in units of GeV), for the branching ratios 



BR mo dci(^fc — ► bb) , 
BR model (h k - W+W-) 



BR mo dcl(^A: — ► T~ + T~ 

BR mo dci(/ifc -> 77) , 



for the following ratios of partonic cross sections 



£model(gg -» h k ) 
gmodel(g^ -» ^fc^' 

q~modei(gg — > ^g) 

a SM (qq ^ HZ) 



O~modcl(bb — » /i fc ) 

0sm(&& -> #) 

gmodel^g ~> ^fcW^" 

o-sM(g'g -> W- 
g G {u,d,c, s,b} , 



O"model(fep — » h k b) 

&SM{bg -> Hb) 
(q,q') E {(u,d), (c,s)} 



and for the hadronic cross section ratio 



Q-modcijpp -> hk via VBF) 
o~sm(pP — » # via VBF) 

for fc 6 {1, . . . ,n H igg S }- 



With this information and some internally provided functions (described be- 
low) HiggsBounds can calculate the same input quantities as required for Input 
Option hadr (subject to some approximations, as described below). It is as- 
sumed that the model cross sections for the processes bg — > hj) and bg — ► h^b 
are equal to each other. 

Hadronic Higgs production cross sections can be written as a sum of cross sec- 
tions of partonic processes a nm -,H+y(S,mH) convoluted by parton luminosity 
functions: 

a(pp -> H + y, m H ) = dr ^ — — & nm ^ H+y (s = tS, m H ) , (41) 

Jto {n,m} dT 

=: ^2 a{pp — > nm —> H + y^mn) ■ 

{n,m} 



Here, the sum includes all distinct parton pairs only once. The centre-of- 
mass energy of the colliding proton (p) and anti-proton (p) is assumed to 
be = 1.96 TeV. The symbol y indicates possible further particles in the 
hard partonic process, like the or Z. The parton luminosity functions are 
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defined as 



AB j. j. 

fn/A{x, fJ-F)fm/B{-, ^f) + fm/A{x, /^F ) /n/B ( ~ , Hf) 



dr Jt x 1 + 8 n 



(42) 



where fn/A(x, Vf) denotes the density of partons of type n in the hadron A 
carrying a fraction x of the hadron's momentum at the scale fip and (A, B) = 
(p,p) for the Tevatron. The lower bound of the r-integration (to) is determined 
by the minimal invariant mass of the parton system, sq = tqS. The ratio 
Qmodei(^) of Eq. (13T|) factorises into a ratio of branching ratio predictions and 
a ratio of hadronic production cross section predictions. The calculation of 
the ratio of branching ratios is straightforward. In the following, we show an 
approximate relation (that in special cases is exact), which allows to calculate 
the ratios of hadronic cross sections for the Higgs production process P, 

Mp) := j (43) 

from the knowledge of ratios of partonic cross sections. Under the assumption 
that the ratio of the partonic cross sections, 

*£»(*,*.„) := iga§^ , (44) 
cr nm ^ H+y {s,m H ) 

has, at most, a mild dependence on the parton-system centre-of-mass energy 
squared, s, we can write Eq. fT4"3l) as 



z?rm~ T}H+y ( - ,<Jsm(pP ^ nm -> H + y,m H ) 
A ) «£„ ( °' H) ^ {p p^H + y, mH ) ' (45) 

with s denoting the partonic production threshold, s = {m H + m y ) 2 , with 
m y = in the case of single Higgs boson production. In this formula, only 
ratios of SM hadronic cross sections are needed, which we provide in our 
program. 

In order to provide the cross section ratios in Eq. (I45p with state-of-the-art 
theoretical accuracy, we pursue the following strategy. The K-factors (NNLO 
QCD + NLO EW) for the WH and ZH production cross sections at the 
Tevatron are, to a very good approximation, independent of the generation of 
the initial state quarks and thus drop out in the ratios. The only non-trivial 
case is the single Higgs production which has contributions from gluon fusion 
and bb annihilation. For the calculation of the gluon fusion cross section, we 
assume a constant NNLO QCD K-factor of 3.6 [51] in the relevant range of 
70 GeV < m H < 300 GeV. In this range, our assumption is consistent with 
the m#-dependent K-factor reported in [51] within the uncertainty estimate 
also given in [5T] . In order to obtain an approximate NNLO QCD cross section 
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Fig. 1. Relative contributions (in %) of different SM partonic processes to cross 
sections of hadronic processes: (a) single Higgs production, (b) HW^ production 
and (c) HZ production. 



prediction for the b b-bar annihilation process, we used the fact that this cross 
section is known to be well approximated by the LO prediction in the MS-bar 
scheme with the factorisation and renormalisation scale choice of |52j . 
The SM cross section ratios from Eq. ( 1451) are then obtained by normalising 
the cross sections for one parton configuration to the appropriate sum of cross 
sections. As a cross-check, we confirmed that our results for Higgs production 
cross sections are in good agreement with the compilation of results by the 
TEV4LHC working group j27j. Fig. [I] shows our result: the relative contribu- 
tions from different parton configurations to the total hadronic cross section 
for single Higgs production, HW production and HZ production in the range 
relevant for the implemented Tevatron analyses. 

Relation (1451) becomes an exact identity for all models which differ from the 
SM only by having different values for coupling constants, i.e. where the ratios 
Rnm V become s-independent. Moreover, for models which give s-dependent 
ratios R*+ y , Eq. ([45]) can still be a good approximation because the luminosity 

i r>pp 

functions are steeply falling functions of r and therefore strongly favour 
the threshold regions of the convolution integrals. However, in this case, the 
user should make sure that Eq. (1451) is indeed a good approximation for the 
model under study before using it. In case of doubt, one can always resort to 
Input Option hadr. 

As an example, we show here how the normalised model cross section Q mo dei(AT) 
is obtained on the basis of Eq. (1451) for the search topology pp — > H — > 
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W + W — » l + l' of Table [2] for a model Higgs boson /i fe : 

^^^ k (so,m hk ) a SM (pp ^bb^ H, m H ) \ ^ 
°bhih k (*o, o"sm(pp -> H,m H ) J 

x BR^f^+iy- (m fefc ) I 



If the model under study predicts kinematical distributions which differ sig- 
nificantly from the distributions of the reference SM processes, the exclusion 
result, returned by HiggsBounds, can only be considered an estimate. Further 
comments regarding this limitation from the end of Section 12.2.11 apply here 
too. 



2.3.3 Tevatron input (HiggsBounds Input Option effC): effective couplings 

The user is asked to provide model predictions for the Higgs boson masses 
m/j fe and total decay widths T tot (hk) (in units of GeV), and for the following 
normalised effective couplings squared: 



' model \ 2 / „modcl\ 2 model\ 2 / -model \ 2 / „model ' 

9h k gg \ I 9h kT r \ I 9h k ZZ \ I 9h k WW \ y h k ff,eS 

y 9ngg J \9h-yjJ \9hzzJ \9hwwJ \ 9 H ff 



(46) 



with k G {1, . . . , ^Higgs} and / G {b, r}. The reference coupling constants gf^ z 

J Hff 



and gf^lf are defined in Eqs. (121)1) and (127j) . respectively, and gfrww 



(9 b H M wwT=[-™w) ■ (47) 



A convenient way to evaluate the ratios of the loop-induced effective couplings 
might be to use the ratio of partial decay widths: 



' „model \ 2 "pmodel /-™ N 
9h k vv' \ _ 1 h k ^vv'\ m h kl 



„SM / PSM l^y, \ 

9Hvv> / 1 H^vv'\" L H I 



(48) 

m H =m hk 



with vv' G {gg, 77}. Of course, the user may apply Eq. (T4"8"l) also for vv' G 
{ZZ, WW} if rriH is large enough or the Ts for off-shell vector bosons are used. 
If the user has only the absolute partial decay widths, T^°^p(mh k ) , available, 
he or she can use the built-in functions for the calculation of SM quantities 
(see Table [12] below for details) in order to calculate the right-hand side of 
Eq. 
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From the input specified in Eqs. (1241) and (|46j) HiggsBounds evaluates all 
partonic cross section ratios, as required in the Tevatron part of Input Option 
part, using the following formulae which are valid in the effective coupling 
approximation: 



&model{bb — > h k ) 



°su{gg 

&modcl(bg 



H) 

h k b) 



a SM (bb^H) & SM (bg ^ Hb) 

CTmodei(qq' ^ h k W + ) _ a mo dci(q'q -> h k W~) 
vsu(qq' HW+) a SM (q'q HW~) 

Vmodei(q"q" -> h k Z) 
vsu(q"q" -> HZ) 

for (q, q') G {(u, d), (c, s)} and g" G {«, c, s, b}. 



' „modcl N 

K 9l^g , 

„model 
gh k bb,cS 
SM_ 
i>Hbb 

( model \ 2 
gh k ww \ 
9hww J 

(„model \ ^ 
<?i§J ' 



(49) 

(50) 

(51) 
(52) 



In order to facilitate the calculation of the ratios of the hadronic cross sections 
for h k production via VBF, we obtained the following SM cross section ratio 
for pp collisions with 1.96 TeV centre-of-mass energy using VBFNLO [29] : 



&sm(pp H via WW fusion) 
°"sm(pp ~~ > H via ZZ fusion) 



3.35 



(53) 



The right-hand side is a mean value. The actual value varies slightly with m# 
(±3.6% in the relevant mass range of 70 GeV < mn < 300 GeV). Using this 
mean and neglecting interference effects, which affect the result below 1% (cf. 
[53]). we can calculate the proportions of WW- and ZZ-fusion contributions 
to the SM VBF Higgs production cross section for pp collisions with 1.96 TeV 
centre-of-mass energy: 



R 



ww . 

VBF 



R 



ZZ 
VBF 



vsmipv ~^ H via WW fusion) 

csm(pp — > H via VBF) 
c"sm(pp — > H via ZZ fusion) 
vsm(pP — > H via VBF) 



77% 



23%. 



(54) 
(55) 



The slight variation of the right-hand side of Eq. ([53 
variation of the above proportions by less than 1%. 



with rriH causes a 



Using the above SM predictions, we can calculate the model predictions for 
the ratios of the hadronic cross sections for h k production via VBF quite 
accurately, as follows: 



o-modeiQpp -» h k via VBF) 
&sm(pP — > H via VBF) 



R 



ww 

VBF 



/ „ model 

/ g>i t ww 



V 9HWW ) 



/ „modcl \ ^ 
, tjZZ / 9h k ZZ \ 



(56) 
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Branching ratios, additional to the ones we have already described when 
discussing the LEP limits in Section I2.2.2[ are evaluated along the lines of 
Eq. (EE): 



BR mo dci(^fc -> 77) 



BRsm(# 



ll){m H ) 



m H =m hk 



X 



X 




(57) 



BR modcl (h k -> WW) 



BRsm(# 



WW)(m H ) 



rtot(^fe) 



m H =m hk 



X 



model 
h k WW 



) 



2 



X 



9hww 



(58) 



At this point, all quantities which are required as Tevatron input by Input 
Option part are evaluated using the effective coupling approximation. The 
evaluation of the model predictions Q model (AT) can then proceed as exemplified 
in Section 12.3.21 above. 

For search topologies, of which the amplitude is proportional to one monomial 
of Higgs coupling constants, effective couplings just rescale all kinematical dis- 
tributions of the reference process by a common factor without changing their 
shape. Consequently, in such cases, limits on topological cross sections can 
be exactly applied. For the application of limits obtained from SM analyses, 
which combine a set of search topologies, the SM-likeness check (see Section 
12. ip ensures that the kinematical distributions arising from models that pass 
this check closely resemble the distributions of the SM. Thus, also in this case, 
limits on topological cross sections can be exactly applied. 

Two slightly less straight-forward cases require further discussion. 

a) For the single Higgs production topology, we add the gluon-fusion and bb 
annihilation contribution together. Firstly, one might expect that the boost 
distribution of the produced Higgs boson is different for the two partonic 
processes and thus the kinematical distributions of the Higgs decay products 
would depend on the relative contributions of the two processes to the signal. 
However, it turns out that the gluon and bottom parton distributions (using 
MSTW 2008 PDFs [SI] ) in the relevant region (lower x values) are, to a very 
good approximation, proportional to each other, i.e. the two parton processes 
do not give rise to different boost distributions of the produced Higgs boson. 
Secondly, it is known that the event selection cuts used in single Higgs analy- 
ses select Higgs bosons with a certain p?, which is known to lead to different 
distributions for processes arising from gluon-gluon and bb initial states. This 
change in the kinematical distribution becomes relevant if the coupling of the 
Higgs to a bottom quark is strongly enhanced such that the bb annihilation 
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competes with gluon-fusion contribution, where the latter dominates in the 
SM. Indeed, it turns out that, for such a scenario, out of the implemented 
set of analyses, the most sensitive are those that involve the search topology 
pp — > H — ► tt [!Mf2"5] . However, these analyses reject events with extra jets of 
large transverse momentum, where the deviation of the j>r-distributions be- 
comes substantial [55]. Therefore, we do not expect a strong influence of the 
moderately changed kinematical distributions on the cross section limits. This 
expectation is supported by the fact that we reproduce the MSSM exclusion 
plots in P^PS] reasonably well. 

b) For Higgs production via vector boson fusion, a deviation from the SM 
proportion of WW and ZZ initiated processes caused by different effective 
couplings g\ zz and g\ ww , does not lead to significant changes in the kine- 
matical distributions compared to the SM reference process. 

In summary, for the implemented set of analyses, we expect the limits on 
topological cross sections to be applicable with high accuracy for models where 
the Higgs sector can be faithfully parametrised with effective couplings. 



3 HiggsBounds Operating Instructions 

There are three formats in which the program HiggsBounds can be used: 

• Library of subroutines 

• Command-line version 

• Online version 

The most widely applicable format of HiggsBounds is the command-line ver- 
sion, since this reads all the model data from text files and thus this model 
data can be generated using any package the user wishes. The library of 
HiggsBounds subroutines allows HiggsBounds to be called within other pro- 
grams. If the user just wishes to check a few parameter points, the online 
version provides quick access to all the functionality of HiggsBounds, without 
the need to install the code. 

The HiggsBounds code, the online version and documentation can all be found 
at the URL www . ippp . dur . ac . uk/HiggsBounds . 

The HiggsBounds code is provided in either Fortran 77 or Fortran 90. Both 
codes provide exactly the same functionality and have exactly the same op- 
erating instructions. In fact, the Fortran 77 and Fortran 90 versions of the 
HiggsBounds subroutines can even be called within codes written in Fortran 
90 and Fortran 77, respectively. Both codes have also been tested with a va- 
riety of Fortran compilers, including the free gnu compilers which accompany 
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most Linux distributions. Therefore, the user may download either code and 
the difference will only be apparent if the user wishes to examine the structure 
of the cod^J. 

The library of subroutines, the command-line version and the online version 
share a common set of features, which we will describe first. We will then give 
operating instructions for each of these three HiggsBounds formats individu- 
ally. 

3.1 Common features: Input 
HiggsBounds requires four types of input: 

• the number of neutral Higgs bosons in the model under study (nH) 

• the set of experimental analyses which should be considered (whichexpt) 

• the theoretical predictions of the model under study (a set of input arrays) 

• the format of these theoretical predictions (whichinput) 

Table H] contains further information on the variable nH, and the possible values 
of whichexpt are described in Table[5j Note that the option whichexpt= ' singH ; 
should only be used if neither processes involving Higgs pair production at LEP 
nor processes involving the hj — *> hihi decay are relevant. However, if these 
conditions are met, this option can save significantly on computing time. 

HiggsBounds expects the theoretical input to be in one of three formats, 
labelled by the variable whichinput. These formats were described in detail 
in Section [2] and are briefly summarised in Table [6) In Table [7] and Table [8] we 
assign names to all of the possible input arrays (each array is defined in terms 
of the notation used in Section [5]). These names will prove useful when we 
describe the input requirements of each version of HiggsBounds individually. 

nH (integer) 

1-9 By default, HiggsBounds can be applied to models with 
between 1 and 9 neutral Higgs bosons. 
(This range can easily be extended by the user if required.) 

Table 4 

The possible values of the variable nH, which labels the number of neutral Higgs 
bosons in the model under study. 



6 The maintenance of two separate codes is primarily intended to provide an effi- 
cient way for the authors to confirm that each update is free from implementation 
errors. 
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whichexpt (character (LEN=5)) 



LandT both LEP and Tevatron analyses 

onlyL only LEP analyses 

onlyT only Tevatron analyses 

singH only analyses for processes involving one Higgs boson 

Table 5 

The possible values of the variable whichexpt, which indicates which subset of 
experimental analyses will be considered by HiggsBounds. 



whichinput (character (LEN=4)) 

effC Masses, total decay widths, 

ratios of effective couplings squared, BR(/ij — ► hihi). 

part Masses, total decay widths, ratios of LEP cross sections, 

mainly ratios of partonic Tevatron cross sections, branching ratios. 

hadr Masses, total decay widths, ratios of LEP cross sections, 

ratios of hadronic Tevatron cross sections, branching ratios. 

Table 6 

The possible values of the variable whichinput, which indicates the format of the 
theoretical predictions provided by the user. (See Section [2] for a more detailed 
description of each of these settings). 
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input arrays (double precision) 

Mh(nH) m hi in GeV 

GammaTotal (nH) r tot (/i,) in GeV 

g 2hjbb(nH) l^ 1 ) , OP = 66 

V 9 H(OP) / 

g2hjtautau(nH) t + t~ 
g2hjWW(nH) WW 
g2hjZZ(nH) ZZ 
g2hjgaga(nH) 77 
g2hjgg(nH) gg 



g2hjhiZ(nH,nH) 



model 
UhjluZ 



BRJijbb(nH) ER(hj — > OP), OP = 66 
BRJijtautau(nH) r + r~ 
BRJijWW(nH) WW 
BRJi j gaga ( nH ) 77 



BRJijhihi(nH,nH) ER(hj -> hihi 



Table 7 

Input arrays for model predictions for effective normalised squared couplings and 
branching ratios recognised by HiggsBounds. The size of each array is given in 
brackets in the first column. See Section [21 for the description of the notation 
used in the second column. The elements of BRJijhihi are ordered such that 
BRJijhihi(j ,i)= BR(hj — >■ /ij/tj). 
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1 Tl TM 1 T Q yy aire C f\ W t 

1I1JJUI aliajo 


(double precision) 




CS_lep_hj Z_ratio (nH) 




P = e+e~ 




fC! "1 or, hnh-i rati nfnH n"H^ 
KjO LcU LI J ±1 J L CLL1U v ■'■■'■■n. j nn J 




e + e~ 


— > h h ■ 

r Ibj llr<l 


CS_tev_pp_hj _ratio (nH) 




PP ~ 


► hi 
3 


tpv r>n hih "rat i n CnH') 

\~i kJ LiC v _UU LI J U L a U 1U V 1J.11 / 




PP ~ 


> hh ■ 

3 


CS_tev_pp_hj W_ratio (nH) 




PP ~ 


> hiW 

3 


tpv nn hiZ "ra1~io(TiH^ 

k_j \j w v i_/ ky l x i i—i j_ u _l w \ xxi x y 




PP 


> hiZ 


CS_tev_pp_vbf _ratio (nH) 




vP 


> /i, viaVBF 

3 


CS_tev_gg_hj _ratio (nH) 




nm = gg 




CS_tev_bb Jij _ratio (nH) 




bb 




tpv nH h i Wn Tat. i o CnH") 

Vy I^J Lr V-/ V LXV-X LX 1 vv L> L CX u _L w \ XXX X / 


r ? h j +W+ 
iMim i 


nm = ud 




vjO _ L, tJ V o LI J W |J L CLU1U V-'-'-f 1 / 




cs 




CS_tev_ud_h.j Wm_ratio (nH) 




nm = du 




CS_tev_cs_h.jWm_ratio (nH) 




sc 




CS_tev_dd_hj Z_ratio (nH) 


„hj+Z 


nm = dd 




CS_tev_mi_hj Z_ratio (nH) 




uu 




CS_tev_ss_hj Z_ratio (nH) 




ss 




CS_tev_cc_hj Z_ratio (nH) 




cc 




CS_tev_bb Jij Z_ratio (nH) 




bb 




CS_tev_bb_hjb_ratio(nH) 




nm = bb 





Table 8 

Input arrays for model predictions for cross section ratios recognised by 
HiggsBounds. The size of each array is given in brackets in the first column. The 
LEP or hadronic Tevatron cross section ratios R a (P) are defined in Eq. &43\) and 
the partonic Tevatron cross section ratios B^m V are defined in Eq. UU). 
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3.2 Common features: Output 



HiggsBounds provides the user with four types of output: 

• whether the parameter point is excluded at the 95% CL. or not (HBresult) 

• the identifying number of the process with the highest statistical sensitivity 
(chan). 

• the number of Higgs bosons which have contributed to the theoretical rate 
for this process (ncombined) 

• the ratio of the theoretical rate Qmodei to the observed limit Q b s for this 
process (obsratio). 

Table M shows the possible values of HBresult and obsratio, which are com- 
plementary. Tables ITTJ1 and [TT1 discuss chan and ncombined respectively. If the 
library of subroutines or the command-line versions are used, the key to the 
process numbers is written in the file Key.dat. In the online version, this 
information appears on the screen. 

HBresult obsratio 

(integer) (double precision) 

> 1.0 parameter point is excluded 

1 < 1.0 parameter point is not excluded 

-1 < 0.0 invalid parameter set 

Table 9 

The possible values of the output variables HBresult and obsratio, which indicate 
whether a parameter point has been excluded at the 95% C.L. by the experimental 
results under consideration. 

chan (integer) 

l-[# of considered analyses] See the file Key.dat for the definition of each 

process number. Key.dat is automatically 
generated when either the command line or the 
subroutine version of HiggsBounds are used. 

Table 10 

Further information about the output variable chan, which stores the reference 
number of the process with the highest statistical sensitivity. 
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ncombined (integer) 



1-iih Number of Higgs bosons which have contributed to the theoretical 

rate for this process. The number depends on delta_Mh_LEP 
or deltaJVIh.TEV. 

Table 11 

Further information about the output variable ncombined. 
3. 3 Library of subroutines 

Installation 

The HiggsBounds code can be compiled to form a library of subroutines using 
the following commands: 

. / configure 
make libHB 

A program which wishes to use the HiggsBounds subroutines can be compiled 
and linked to the library by adding -L<HBpath> -1HB to the command line, 
for example, 

gfortran myprog.f90 -o myprog -L<HBpath> -1HB 
where <HBpath> is the location of the HiggsBounds library. 

Subroutine initialize_HiggsBounds 

The subroutine initialize_HiggsBounds must be called before any other 
HiggsBounds subroutine. It performs some preparatory operations such as 
reading in the tables of data. It is called as: 

call initialize_HiggsBounds(nH, whichexpt) 

When using the HiggsBounds subroutines in another code, the subroutine 
initial ize_HiggsBounds must be called only once, before any other 
HiggsBounds subroutine is called. 

Subroutines run_HiggsBounds_ef f C , run_HiggsBounds_part and 
run_HiggsBounds_hadr 

The subroutines run_HiggsBounds_ef f C, run_HiggsBounds_part and 
run_HiggsBounds_hadr perform the main part of the HiggsBounds calcula- 
tions. They set the value of whichinput to be ef fC, part and hadr respec- 
tively and therefore require different arguments. These subroutines are called 
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as: 



call run_HiggsBounds_ef f C (nH , Mh , GammaTotal , 

& g2h j bb , g2h j t aut au , g2h j WW , g2h j ZZ , 

& g2hjgaga,g2hjgg,g2hjhiZ, 

& BRJijhihi, 

& HBresult , chan, 

& obsratio, ncombined ) 

call run_HiggsBounds_part (nH,Mh, 

& CS_lep_hjZ_ratio , CS_lep_hjhi_ratio , 

& CS_tev_gg_hj_ratio , CS_tev_bb_hj_ratio , 

& CS_tev_bb_hjb_ratio , 

& CS_tev_ud_hj Wp_ratio , CS_tev_cs_hj Wp_ratio , 

& CS_tev_ud_hj Wm_ratio , CS_tev_cs_hj Wm_ratio , 

& CS_tev_dd_hjZ_ratio , CS_tev_uu_hjZ_ratio , 

& CS_tev_ss_hjZ_ratio , CS_tev_cc_hjZ_ratio , 

& CS_tev_bb_hjZ_ratio , 

& CS_tev_pp_vbf _ratio , 

& BR_hjbb,BR_hjtautau, 

& BR_hjWW,BR_hjgaga, 

& BR_h.jh.ihi, 

& HBresult , chan, 

& obsratio, ncombined ) 

call run_HiggsBounds_hadr (nH,Mh, 

& CS_lep_hjZ_ratio , CS_lep_hjhi_ratio , 

& CS_tev_pp_hj_ratio, CS_tev_pp_hjb_ratio , 

& CS_tev_pp_hjW_ratio , CS_tev_pp_hjZ_ratio , 

& CS_tev_pp_vbf _ratio , 

& BR_hjbb,BR_hjtautau, 

& BR_hjWW,BR_hjgaga, 

& BRJijnihi, 

& HBresult , chan, 

& obsratio, ncombined ) 



Each of these arguments must be supplied. However, if a branching ratio, 
effective coupling or cross section is believed to be irrelevant, the corresponding 
array may be filled with zeros. This will ensure that the value of Q mo dei f° r 
processes involving this quantity will also be zero. For example, in the MSSM, 
processes involving the decay hj — > 77 will rarely be the process with the 
highest statistical sensitivity of Qmodci/Qcxpec and, consequently, it may be 
convenient to set the arrays g2hjgaga and BRJijgaga to zero for simplicity. 
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Also, depending on the value given for whichexpt, some of the input arrays 
will be ignored within HiggsBounds. For example, if whichexpt= < onlyT , , the 
branching ratio for the Higgs cascade decay hj — > hihi will not be relevant. 
Therefore, setting this array to zero in this case will not affect the HiggsBounds 
results in this case. 

Subroutine f inish_HiggsBounds 

The subroutine f inish_HiggsBounds should be called once at the end of the 
program, after all other HiggsBounds subroutines^! It is called as: 

call f inish_HiggsBounds 

Functions for Standard Model branching ratios, total decay width and cross 
sections 

The HiggsBounds library also allows users access to the Standard Model Higgs 
branching ratios, total decay width and production cross sections, which are 
used internally by HiggsBounds. We use Standard Model Higgs branching 
ratios and total decay width from the program HDECAY 3.303 [26]. The SM 
hadronic cross sections have been obtained from the TEV4LHC Higgs Working 
Group [27] (see Table [T2l for references to the original works) with the exception 
of the a SM (pp — > bg — > bH) cross section. The latter cross section has been 
calculated with the program HJET 1.1 [28] for a set of different cuts on the 
transverse momentum and pseudo-rapidity of the fe-quark, which are needed 
internally in order to apply correctly the results of some Tevatron analyses 
and which were not available from [27]. From this set, only the cross section 
without cuts is externally provided. 

Examples 

We have provided three example programs which demonstrate the use of the 
HiggsBounds subroutines. The first example relates to the Fourth Genera- 
tion Model and is contained in the file example-SM_vs_4thGen.F (see Fig. H] 
below for a code listing). This program uses the HiggsBounds functions for 
the SM branching ratios and SM total decay width to calculate the Higgs 
decay width and the effective normalised squared couplings in the SM and a 
simple Fourth Generation Model. This information is then used as input for 
the subroutine run_HiggsBounds_ef f C, which is called once with SM input 
and once with Fourth Generation Model input. The example will be described 
in more detail in Section 14.21 Once the HiggsBounds library has been com- 
piled (using ./configure ; make libHB as described previously), the code 



7 In the Fortran 90 version of the code, the subroutine f inish_HiggsBounds is used 
to deallocate the allocatable arrays used within HiggsBounds. 



37 



function 


(double precision) 








SMGammaJi(Mh) 


rf5f(fc) 






126] 


SMBR_Hbb(Mh) 


BR SM (H -> OP), OP= 


bb 




[26] 


SMBR_Htautau(Mh) 




T~T + 




I2S] 


SMBR_HWW(Mh) 




WW 




[26] 


SMBR_Hgamgam(Mh) 




77 




[26] 


SMBR_Hgg(Mh) 




fjfj 




[26] 


SMBR_HZgam (Mh) 




Z 1 




[26] 


SMBR_HZZ(Mh) 




zz 




[26] 


SMBR_Htoptop (Mh) 




tt 




[26] 


SMBR_Hcc(Mh) 




ce- 




[26] 


SMBR_Hss (Mh) 




ss 




[26] 


SMBR_Hmumu (Mh) 




(jTn + 




[26] 


SMCS_tev_pp_qq_HW (Mh) 


a SM (P), P= 


PP -» QQ 


-» -H"W 


[56] 


SMCS_tev_pp_qq_HZ (Mh) 




PP ~ *• QQ 


-> 


[56] 


SMCS_tev_pp_gg_H (Mh) 




pp^ gg 


iJ 


|50|5l] 


SMCS_tev_pp_bb_H (Mh) 




pp — > bb 


-► if 


[52] 


SMCS_tev_pp_vbf _H (Mh) 




pp ^ H 


via VBF 


[57] 


SMCS_tev_pp_bg_Hb (Mh) 




pp -> bg 


-> 


[28] 



Table 12 

Standard Model branching ratios, total decay widths in units of GeV and hadronic 
Tevatron cross sections in units of pikobarn provided as functions by HiggsBounds, 
together with references. Each function takes a Higgs mass Mh (double precision) 
as its argument. 



example-SM_vs_4thGen.F can be compiled and run with the commands: 

gfortran example-SM_vs_4thGen.F -o HBfourthgen -L<HBpath> -1HB 
. /HBfourthgen 

where <HBpath> is the location of the HiggsBounds library. 

The files HBwithFH.F and HBwithCPsuperH.f demonstrate the use of 
HiggsBounds subroutines with the publically available programs FeynHiggs [T3] 
and CPsuperH [18], respectively. We refer the reader to the extensive comments 
contained within these example files for further details. 
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3.4 Command line version 
Installation 

In order to be able to call HiggsBounds from the command line, it should be 
compiled using the commands 

. / configure 
make 

Command line and input file format 

In the command-line usage of HiggsBounds, the arrays containing the theoret- 
ical model predictions are read from text files. The other options are specified 
in the command line, which is of the form: 

./HiggsBounds <whichexpt> <whichinput> <nH> <prefix> 

The variable <pref ix> is a string which is added to the front of input and 
output file names and may include directory names or other identifying infor- 
mation. 

Table [13] describes the contents of each input file. Note that each input file 
should start with a line number. The input files should not contain any com- 
ments or blank lines. The line number identifies the predictions which belong 
to the same model parameter point in different files. 

Care should be taken with the order of the array elements in the files. This is 
best illustrated by an example, where we will use tih — 3. The one dimensional 
arrays, e.g. Mh, should be given in the order 

Mh(l), Mh(2), Mh(3) 

However, not all of the elements of the two dimensional arrays are required. 
Only the lower left triangle (including the diagonal) is required from the arrays 
g2hjhiZ and lepCS_hjhi_ratio, since they are symmetric, e.g. 



g2hjhiZ(l,l) g 2hjhiZ(l,2) g 2hjhiZ(l,3) 
g 2hjhiZ(2, 1) g 2hjhiZ(2,2) g 2hjhiZ(2,3) 




i.e. the elements in the input file should be written in the order 



g2hjhiZ(l,l), g2hjhiZ(2,l), g2hjhiZ(2,2), g 2hjhiZ(3, 1), 
g2hjhiZ(3,2), g2hjhiZ(3,3) . 
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file name data format 



MH_GammaTot . dat 




k, 


Mh, GammaTotal 


effC.dat 




k. 


g2h j bb , g2h j t aut au , g2h j WW , g2h j ZZ , 








g2rijgaga,g2hjgg, 








some elements of g2hjhiZ 








(lower left triangle - see example) 


LEP_HZ_CS_ratios . dat 




k. 


CS_lep_hj Z_ratio 


LEP_2H_CS_ratios . dat 






some elements of CS_lep_hjhi_ratio 








(lower left triangle - see example) 


TEV_H_Ojet_partCS_ratios . 


dat 


k, 


CS_tev_gg_hj _ratio , CS_tev_bb_hj _ratio 


1 J_j V _T1_1 J fci L_pd± LUo_i dL-lOo . 


Cld. L 


k CS_tev_bb_hjb_ratio 


TEV_HW_partCS_ratios . dat 




k. 


CS_tev_ud_hj Wp_rat io, 








CS_tev_cs_hjWp_ratio, 








CS_t e v_ud_h j Wm_r at i o , 








CS_tev_cs_hjWm_ratio 


TEV_HZ_partCS_ratios . dat 




k. 


CS_tev_dd_hjZ_ratio, 








CS_tev_uu_hjZ_ratio , 








CS_tev_ss_hjZ_ratio, 








CS_tev_cc_hjZ_ratio , 








CS_tev_bb_hjZ_ratio 


TEV_H_vbf_hadCS_ratios . dat 


k, 


CS_tev_pp_vbf _ratio 


TEV_lH_hadCS_ratios . dat 




k, 


CS_tev_pp_hj_ratio, 








CS_tev_pp_hjb_ratio, 








CS_tev_pp_hjW_ratio, 








CS_tev_pp_hjZ_ratio , 








CS_tev_pp_vbf _ratio 


BR_1H . dat 




k. 


BRJij bb , BRJi j t aut au , BRJi j WW , BRJi j gaga 


BR_2H . dat 




k, 


some elements of BR_h.jh.ihi 








(row by row, without diagonal 








- see example) 


additional . dat (optional) 




k, 





Table 13 

Names and data format of all HiggsBounds input files. The right column shows 
the order of the input data arrays within one line of the input file. For the or- 
der within the arrays, see example, k is the line number. Note that the array 
CS_tev_pp_vbf _ratio appears in two different input files. However, these files will 
never be required by HiggsBounds simultaneously. 
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For the array BRJijhihi, only the off-diagonal components are required 

I3RJijhihi(l, 1) BRJijhihi(l, 2) BRJijhihi(l, 3)^ 
BRJijhihi(2, 1) BR_hjhihi(2, 2) BRJijhihi(2, 3) 
J3RJijhihi(3, 1) BR_hjhihi(3, 2) BRJijhihi(3, 3) 

since the diagonal elements are not physical quantities. Therefore, the elements 
should be written in the order 

BRJijhihi (1,2), BRJijhihi (1,3), BRJijhihi ( 2, 1), BRJijhihi (2, 3), 
BRJijhihi (3,1), BRJijhihi (3, 2) 

in the input file. 

The file additional . dat is optional. If it is included, it can have any number 
of columns greater than 1 (as for the previous files, the first entry on each 
line should be the line number). It is envisaged that this input file will be 
particularly useful when parameter scans are performed over a variable which 
is not required by HiggsBounds but helpful plotting the results. For example, 
in the case of the complex MSSM, additional . dat could be used to store the 
values of tan (5 and the charged Higgs mass. 

As in the subroutine version, the command line version of HiggsBounds ex- 
pects a subset of the total list of input arrays, which depends on the cho- 
sen setting of whichinput. The maximal list of files used for each value of 
whichinput is given in Table [5) 

As discussed for the subroutine version, some of the arrays will not be rel- 
evant for some of the choices for whichexpt. The command line version of 
HiggsBounds will consider the list of input files appropriate to the setting 
whichinput and then only attempt to read any of these input files if the value 
chosen for whichexpt means that at least one of the arrays it contains will 
be directly used. Table [15] contains a list of which input files are actually rele- 
vant to each value of whichexpt. For example, if whichinput = 'hadr', and 
whichexpt = 'LandT', then HiggsBounds requires the input files: 

MH_GammaTot . dat, BR_2H . dat, BR_1H . dat, LEPJlZ_CS_ratios . dat, 
LEP_2H_CS_ratios . dat, TEV_lHJiadCS_ratios . dat . 

However, if whichinput = 'hadr' and whichexpt = 'onlyL', HiggsBounds 

requires the input files: 

MH_GammaTot . dat, BR_2H . dat, BR_1H . dat, LEPJlZ_CS_ratios . dat, 
LEP_2H_CS_ratios.dat . 
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whichinput = 'part', 



hadr ' , 



effC 



MH_GammaTot . dat 



MH_GammaTot . dat 



MH_GammaTot . dat 



BR_2H . dat 



BR_2H . dat 



effC.dat 



BR_1H . dat 



BR_1H . dat 



BR_2H . dat 



LEP_HZ_CS_ratios . dat 



LEP_HZ_CS_ratios . dat 



additional . dat 



LEP_2H_CS_ratios . dat 



LEP_2H_CS_ratios . dat 



TEV_H_0 j et_partCS_ratios . dat TEV_lHJiadCS_ratios . dat 
TEV_H_1 j et_partCS_ratios . dat additional . dat 
TEV_HW_partCS_ratios . dat 
TEV_HZ_partCS jratios . dat 
TEV_H_vbf _hadCS_ratios . dat 
additional . dat 



The list of possible input files for each value of whichinput. Note that some input 
files may not be relevant, depending on the value of whichexpt. In this case, they 
are not required. See Table\T5\for more details. 



As a third example, if whichinput = 'hadr' and whichexpt = 'onlyT', 
HiggsBounds requires the input files: 

MH_GammaTot.dat, BR_lH.dat, TEV_lH_hadCS_ratios.dat . 

In each of these three examples, HiggsBounds will also read the file 

additional .dat 
if it exists. 

As for the subroutine version, if the user does not require processes involving a 
particular branching ratio or cross section ratio to be checked by HiggsBounds, 
that particular array can be filled with zeros. 

Output file format 

When the command line version of HiggsBounds is used, the output is writ- 
ten to the file <pref ix>HiggsBounds_results .dat. A sample of the out- 
put is shown in Fig. [2j The key to the process numbering is written to 
<pref ix>Key . dat. 



Table 14 
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name of input file values of whichexpt 

which this file is relevant to 





T anrIT 


onlyL 


onlyT 


o i n 1 — 1 

Slllg-tl 


MU P -i mm -iTa+- A -3 +■ 


y 


y 


y 


y 


eilL' . CL3.X 


y 


y 


y 


y 


T T7"D U7 f"*Q vtH — i /~\ f-i 


y 


y 




y 


LhF_zn_Cb_ratios . dat 


y 


y 






lhVJi_Uj et_partCb_ratios . dat 


y 




y 


y 


lEjV_rl_l J GX_pajrXL»o_raT,10S . 0.3X 


y 




y 


y 


ihv_HW_partCb_ratios . dat 


y 




y 


y 


TEV_HZ_partCS_ratios . dat 


y 




y 


y 


TEV_H_vbf _hadCS_r at i o s . dat 


y 




y 


y 


TEV_lH_hadCS_ratios . dat 


y 




y 


y 


BFL1H . dat 


y 


y 


y 


y 


BR_2H . dat 


y 


y 






additional.dat (optional) 


y 


y 


y 


y 



Table 15 

List of input files, specifying which values of whichexpt each input file is relevant 
to (marked by 'y'). 

Examples 

The HiggsBounds package includes a full set of sample input files for the case 
•an — 3, contained in the folder example. Each filename is prefixed with nH3_. 
To run the command-line version of HiggsBounds with these files as input, 
use, for example, 

./configure 
make 

./HiggsBounds LandT effC 3 ' example/nH3_ ' 

where the values of whichexpt and whichinput can be varied as desired. 
3.5 Website version 

The website allows the user to select the required number of neutral Higgs 
bosons and then generates a html form accordingly. The values of whichinput 
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# generated with HiggsBounds on 31.10.2008 at 11:18 

# settings: LandT, effC 



# column abbreviations 



Mh(i) 

HBresult 

chan 

obsratio 

ncomb 

additional 



line id of input 
Higgs boson masses 

scenario allowed flag (1: allowed, 0: excluded, -1: unphysical) 

most sensitive channel (see below) . chan=0 if no channel applies 

ratio [sig x BR] _model/ [sig x BR]_limit (<1: allowed, >1: excluded) 

number of Higgs bosons combined in most sensitive channel 

optional additional data stored in <pref ix>additional .dat (e.g. tan beta) 



channel numbers used in this file 



3 
4 

124 
134 
157 



(ee)->(h3)Z->(b b)Z (LEP table 14b) 
(ee)->(hl)Z->(tau tau)Z (LEP table 14c) 
(pp)->W(hl)->l nu (b b) (CDF Note 9463) 

(pp) ->h2->tau tau (\protect\vrule widthOpt\protect\href {http : //arxiv . org/abs/0805 . 2491}{arXiv : 0805 . 2491} 
(pp)->hl+... where hi is SM-like (\protect\vrule widthOpt\protect\href {http : //arxiv. org/abs/0804 . 3423}{ar 



(for full list of processes, see Key. dat) 



#cols : 
# 



n 


Mh(l) 


Mh(2) 


Mh(3) 


HBresult 


chan 


obsratio 




ncomb 


additional(l) 


1 


359.121 


271.963 


134.929 


1 


134 


0.212206E- 


-03 


1 


0.246862 


2 


75.0123 


92.8677 


71.9716 


1 


4 


0.306172E- 


-01 


1 


0.714964 


3 


136.293 


345.483 


330.026 


1 


124 


0.640713E- 


-01 


1 


0.434594 


4 


111.377 


220.765 


51.7469 


1 


3 


0.162811 




1 


0.727173 


5 


186.131 


355.002 


146.448 





157 


15.2354 




1 


0.230522 



Fig. 2. Sample output file (written to <pref ix>HiggsBounds_Results . dat ) 

and whichexpt can be chosen and the appropriate theoretical input entered. 
HiggsBounds will then be called with these settings and the result outputted 
to screen. The website version contains the additional feature that it notifies 
the user about the processes with the second and third highest statistical 
sensitivities and the values of obsratio for these processes. This is designed 
to give guidance to the user who, for example, wishes to find an excluded 
region iteratively by adjusting the input quantities. 

The website also contains a selection of pre-filled html forms as examples, 
including entries for the Standard Model, the fermiophobic Higgs model and 
the MSSM with real and complex parameters. 



4 Examples of use 



4-1 Re- evaluation of SM exclusion limit with an improved SM prediction 



Recently, an improved prediction for the Higgs-boson production via gluon 
fusion appeared [58], which incorporates previously unknown mixed QCD- 
electroweak radiative corrections to this process as well as other known rele- 
vant corrections [59J. It turns out that those corrections are relevant, increasing 
the SM cross section prediction by 7-10% compared to the cross section pre- 
diction currently used by the Tevatron collaborations, when the MRST 2006 
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125 130 135 140 145 150 155 160 165 170 175 180 185 190 195 200 125 130 135 140 145 150 155 160 165 170 175 180 185 190 195 200 

m„ [GeV] m„ [GeV] 

Fig. 3. Ratio Qobs/QsM for the most sensitive channel as a function of the Higgs 
mass uih'- (a) The current best available SM combined analyses of CDF and D0 
(solid line) is compared to a re-evaluation using the recently improved SM prediction 
for the gluon fusion process with MRST 2006 PDFs (dashed line) and MSTW 2008 
PDFs (dotted line) \58\ : (b) The impact of the PDF uncertainty on the exclusion 
limit using the improved SM prediction J58| is shown. 

NNLO parton distribution functions (PDFs) [60] are used. However, very re- 
cently, an updated set of PDFs, MSTW 2008 [51], appeared which show a 
decrease in the low-x gluon distribution. In an updated version of Ref. [58] 
the combined effect of using the improved prediction for the SM cross section 
and the latest set of MSTW 2008 NNLO PDFs has been found to result in 
a downward shift of the SM cross section prediction in the region where the 
Tevatron searches have the highest sensitivity by 4-6% with respect to the 
cross section prediction used by the Tevatron collaborations. 

With HiggsBounds the impact of such changes in the SM cross section predic- 
tion on the SM exclusion limit can easily be assessed. Fig. [3] shows the result 
of a re-evaluation of the Tevatron exclusion limits. In Fig. |3^a), the limit 
obtained using the cross section predictions of [5H] with MRST 2006 NNLO 
PDFs and MSTW 2008 NNLO PDFs are compared with the published lim- 
its from the CDF and DO SM combined analyses ([H] for m# < 155 GeV 
and [49j for m# > 155 GeV) which do not include this improved prediction. 
While in the published Tevatron result the cross section limit just touches the 
SM prediction for m H ^170 GeV, the updated SM cross section prediction 
exceeds the cross section limit within an interval of about 5 GeV if MRST 
2006 PDFs are used. However, when the most up-to-date MSTW 2008 PDFs 
are used, the cross section prediction is lower than the one obtained by the 
Tevatron collaborations for m# above about 160 GeV, resulting in a reduc- 
tion of the exclusion power. With the updated prediction, no value of the SM 
Higgs boson mass is excluded at the 95% C.L. Fig. E](b) shows how the com- 
parison between the SM cross section prediction and the exclusion limits from 
the Tevatron is affected by PDF uncertainties. The solid line shows the result 
based on the updated SM cross section prediction using MSTW 2008 PDFs, 
while the dashed lines indicate the PDF uncertainties as given in Ref. |58j . 
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It can be clearly seen that currently no firm exclusion of any SM Higgs mass 
values at the 95% CL. can be established from the exclusion bounds obtained 
at the Tevatron once PDF uncertainties on the SM cross section prediction 
are taken into account. 

4-2 Fourth Generation Model 

A very simple example of physics beyond the SM is a model which extends the 
SM by a fourth generation of heavy fermions [01]. In particular, the masses of 
the 4th generation quarks and leptons are assumed to be (much) heavier than 
the mass of the top-quark. In this case, the effective coupling of the Higgs 
boson to two gluons is three times larger than in the SM. No other coupling, 
relevant to LEP and Tevatron searches, changes significant lj@. Essentially, 
only the partial decay width T(H — > gg) changes by a factor of 9 and, with 
it, the total Higgs width and therefore all the decay branching ratios [52] . 

We can easily calculate the new total decay width and the relevant decay 
branching ratios using the above information and the built-in functions to 
calculate SM quantities using the relations: 

T SM (H - gg) = BR SM (H - gg) T^(H) , 
r m0 dci(# -> gg) = 9 T SM (H -> gg) , 

r t m ? dcl (# ) = T™(H) - T SM (H -> gg) + r model (# - gg) . 

It is then very simple to test the model with HiggsBounds using the effective 
coupling input. 

In Fig. H] we show a Fortran 77 code example which uses HiggsBounds for 
the calculation of the Higgs mass dependence of the ratio Qmodci/Qobs for 
the SM and our simple Fourth Generation Model scenario. From the results 
returned by the subroutine run_HiggsBounds_ef f C, the set (mh, HBresult, 
chan, obsratio) is stored in a file, separately for the two models. 

The results obtained by this sample program are presented in Fig. El which 
shows the ratio (Qmodei/Qobs)" 1 = [cr x BR] obs . ii mit /[cr x BR] modcl as a function 
of the Higgs mass m# for the SM (solid lines) and the Fourth Generation 
Model (dashed lines). For the SM cross section we use here the prediction 
employed by the Tevatron collaborations, see the discussion above concerning 

8 The effective Higgs-photon-photon coupling is suppressed by a factor of about 
(H)/T^(H) for a light Higgs H [62]. The exact suppression factor is irrel- 
evant here, as the Higgs decay into 77 can only lead to relevant search topologies 
at the Tevatron if it is enhanced by more than an order of magnitude compared to 
the SM. 



16 



program SM_vs_f ourth_generation_model 
implicit none 

integer nH , HBresult , chan , ncombined 

real*8 obsrat io , SMGammaTotal , GammaTotal , SMgammaJigg , 
& GammaJigg , SMGammaJi , SMBR_Hgg , Mh , g2h j bb , g2h j t aut au , 
& g2h j WW , g2h j ZZ , g2h j gaga , g2h j gg , g2h j hiZ , BRJi j hihi 

nH=l 

call initialize_HiggsBounds (nH, 'LandT') 
open ( 10, f ile =) example-SM-results . dat ' ) 
open (20 , f ile= ' example-4thGen-results . dat ' ) 

do Mh=90,200,5 
SMGammaTotal=SMGammaJi(Mh) 
SMgamma_hgg=SMBR_Hgg(Mh) *SMGammaTotal 
Gamma _hgg=9dO*SMBR_Hgg (Mh) *SMGammaTotal 
GammaTotal=SMGammaTotal-SMgamma_hgg+GammaJigg 

g2hjbb=ld0 

g2hjtautau=ld0 

g2hjWW=ldO 

g2hjZZ=ldO 

g2hjgaga=ld0 

g 2hjgg=ld0 

g2hjhiZ=0d0 

BR_hjhihi=0d0 

call run_HiggsBounds_ef f C (nH , Mh , SMGammaTotal , 
& g2h j bb , g2h j t aut au , g2h j WW , g2h j ZZ , g2h j gaga , g2h j gg , g2h j hiZ , 
& BRJij hihi , HBresult , chan, obsratio , ncombined) 

wr ite ( 10 , *) Mh , HBresult , chan , obsrat io 

g2hjgaga=ldO/9dO*GammaTotal/SMGammaTotal 
g 2hjgg=9d0 

call run_HiggsBounds_ef f C (nH , Mh , GammaTotal , 
& g2h j bb , g2h j t autau , g2h j WW , g2h j ZZ , g2h j gaga , g2h j gg , g2h j hiZ , 
& BRJij hihi , HBresult , chan, obsratio , ncombined) 
write (20 , *) Mh , HBresult , chan , obsrat io 
enddo 

close(lO) 
close(20) 

call f inish_HiggsBounds 
end 

Fig. 4. Fortran77 code example-SM_vs_4thGen.F which calculates the ratio 
Q model/ Qobs (obsratio) for the most sensitive Higgs search topology for the SM 
and a simple Fourth Generation Model as a function of the Higgs mass (Mh). This 
program has been used to produce the plot shown in Fig. 
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m H [GeV] 

Fig. 5. Ratio Qobs/Q model for the most sensitive channel as a function of the Higgs 
mass ran '■ 4th Generation Model versus SM. The colours indicate whether the most 
sensitive search channel is from LEP (lighter grey) or the Tevatron (darker grey). 

an update of this prediction. Exclusion at 95% C.L. is given if the value of 
(Qmodei/Qobs) -1 drops below 1. Both curves have sections in red colour (lighter 
grey), where the LEP results provide the most sensitive search topology, and in 
blue colour (darker grey), where the most sensitive search topology is provided 
by the Tevatron results. As expected, the LEP results have the highest sensi- 
tivity in the low-mass region up to the kinematical limit of about 115 GeV, 
while for higher mass values the Tevatron results take over. 

The SM results show the familiar exclusion of Higgs masses below « 115 GeV 
from LEP and the small excluded region at about m# = 170 GeV from the 
latest Tevatron SM combined analysis [49] using the results of the slightly 
more conservative results obtained by the Bayesian method (see the discussion 
above for an update of the cross-section prediction employed here). For the 
Fourth Generation Model the LEP exclusion is a bit weaker than for the SM, 
because at low Higgs mass the larger branching ratio for Higgs decays into 
gluon pairs diminishes the bb and r + r~ branching ratios. In the region above 
115 GeV, the most sensitive search topology for the Fourth Generation Model 
Higgs is pp — > H — > W + W~ — > throughout the displayed range. Because 
of the enhancement of the signal cross section due to the enhanced coupling of 
the Higgs boson to gluons, the region excluded by the Tevatron results extends 
from about 145 GeV to the end of the displayed range. Because of the SM- 
likeness criterion, the single topology analyses [39E0] apply, instead of the 
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SM combined analysis, for the Fourth Generation Model. Disregarding this 
criterion would change the (Q m odei/Qobs) _1 values in the Tevatron exclusion 
region to almost exactly 1/9 times the SM values. 



4.3 MSSM 
Theoretical background 

The MSSM Higgs sector consists of two scalar weak isospin doublets and 
possesses therefore five physical Higgs bosons. At tree-level these are the CP- 
even h and H (with irih < m#), the CP-odd A and the charged H^. The 
Higgs sector of the MSSM can be expressed at lowest order in terms of Mz, 
Ma (or M H ±) and tan/3 = v 2 /vi, the ratio of the vacuum expectation values 
of the two doublets. All other masses and mixing angles can therefore be 
predicted. The tree- level masses and couplings are strongly affected by higher- 
order corrections [63]. In the real MSSM (rMSSM), i.e. the model without 
CP-violation, the higher-order corrected masses are denoted as Mh, Mh and 
Ma- In the MSSM with complex phases (cMSSM) the three neutral Higgs 
bosons can mix and the masses are denoted as M hl , M h2 , M h3 (with M hl < 
M h2 <M hs ). 

In order to ensure that the external Higgs bosons have the correct on-shell 
properties, the tree-level Higgs states have to be supplemented with Higgs- 
propagator corrections (which are often the dominant corrections). This is 
done with the help of the Z matrix [TTfM] . a non-unitary matrix, which mixes 
between the Higgs tree-level mass eigenstates (in the order h, H, A) and the 
loop-corrected mass eigenstates 







M 


h 2 


= Z 


H 









In the rMSSM, setting hi = h, h 2 = H, h% = A, one has 

Z31 = Z32 = Z\z = Z23 = 0, Z33 = 1. (60) 

In the following we give the formulas valid for the more general cMSSM, where 
the special case of the rMSSM can be derived with the help of Eq. ( !60|) . The 
normalised cross sections o~f t °^ m , <r^°™, defined through Eqs. (fTj), (jig) and (|20|) 
as part of the input needed by HiggsBounds, can then be found using the 
elements of the matrix Z 



49 



<T = \Zn sm(/5 - a) + Z i2 cos(/3 - a) | 2 

<X = cos(/5 - a) - Z j2 sin(/3 - a)) (61) 

-Zjs(Zn cos(/3 - a) - Z i2 sin(/5 - a))| 2 , (62) 

where a is the angle that diagonalises the CP-even Higgs sector at tree-level. 
It should be noted that for two identical Higgs bosons this leads to cr£°™ = 
and that the following approximate relations hold: 



+ + (63) 

(64) 

where i,j,k are all different!" 9 "]. 



Limits for the real MSSM 

In Ref. [12] the final results for the LEP Higgs searches in the MSSM were 
published. In the case of the rMSSM they had been obtained in four bench- 
mark scenarios: the m™ ax , no-mixing, small a e ff and gluophobic Higgs sce- 
narios (including also some variations), see Ref. [65] for a motivation and 
detailed definition of these scenarios. The computer code used at that time 
was FeynHiggs 2 . 0, and in most cases the top-quark mass was fixed to m t = 
174.3 GeV. The limits had been obtained by a combination of all available 
channels. 

Here we present the limits obtained by HiggsBounds using the latest version 
(2.6.4) of FeynHiggs to calculate the required model input and choosing dif- 
ferent values for m t . Note that the limits given here are in the narrow- width 
approximation, as discussed in Section 12.11 

In the left column of Fig. [6] we show the LEP and Tevatron exclusion bounds 
as obtained with HiggsBounds in the m™ ax (a), no- mixing (b), small a c s (c) 
and gluophobic Higgs scenario (d) in the M^-tan f3 plane for the current top- 
quark mass value, m t = 172 .4 GeV [66J. In the right column of Fig. [6] we show 
the most sensitive channels for each point in the M^-tan j3 plane for the same 
four scenarios (e)-(h) (the colour code is given in the caption). 

The four scenarios show similar features in the excluded regions, as presented 
in the left-hand column. Very low Ma values are excluded by LEP searches 
for all tan/3 values. This is due to the channel e + e~ — > HA — > bbbb as can 
be seen in the right column of Fig. [6j In this region the hZZ coupling is 
suppressed. Only at very low tan (3 values the channel e + e~ — > hZ, h — > bb has 



9 Note that these relations become equalities if a unitary Higgs mixing matrix is 
used. See, e.g., Ref. [T7] . 
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Fig. 6. Left column [(a),(b),(c),(d)]: 95% C.L. level excluded regions in the Ma -tan (3 
plane for the the four rMSSM benchmark scenarios. The colours indicate exclu- 
sion from LEP (M) or the Tevatron (M). Right column [(e),(f),(g),(h)]: chan- 
nel of highest statistical sensitivity: yellow ( ): e + e~ — > hZ,h — > bb; cyan (M) 
e+e" -> HZ,H -* bb; cyan (B): e+e" -> hZ/HZ,h/H -> bb; orange (M/B) 
e+e" -> hA -> bbbb/bbr+T-; green (M): pp -> h/A t+t~ J2§f; blue (U) 
pp ->• h/H/A -> r + r" 125$; magenta (M): pp -> if/A ->• r + r" {2^; magenta (□) 
-> A -> r+r" {2>f; red (■ j: pp -»• W -> 6Wi/ J2>f; Wack (■ ): pp -> if /A r+r" 
J24f ; blue (13): pp — > ii — > W + W~ )39| . Tie contour borders the excluded regions. 
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the highest statistical sensitivity, but still with a somewhat suppressed hZZ 
coupling. At around Ma ~ 100 GeV the kinematical limit for the HA channel is 
reached and other channels take over. There the exclusion drops to a relatively 
low excluded tan (3 value, slightly different in the four scenarios. At around 
Ma ~ 115 GeV the exclusion in the e + e" — ► hZ, h — > bb channel becomes SM- 
like, and Higgs bosons up to the kinematical limit can be excluded by LEP 
searches, reaching tan/3 = 8, 30, 10, 50 in the m™ ax , no-mixing, small a c s and 
gluophobic Higgs scenario, respectively. In the latter one also larger values 
could be excluded, but the Tevatron channels become more important, see 
below. For larger Ma values the LEP exclusion drops to a lower tan/5 value 
between 3 and 11, depending on the scenario. These tan/3 values correspond 
to Mh ~ 208 GeV — Mz, i.e. the kinematical limit and are excluded for all 
Ma values above ~ 115 GeV. 

The Tevatron searches exclude a parameter space, depending on the scenario, 
between Ma = 110 GeV and 170 GeV, and between tan/3 = 39 and 60 (where 
we stopped our scan, but in principle larger values can be excluded). The 
relevant channels are pp — > h/H/A + X, h/H/A — > t + t~, see the right-hand 
column of Fig. The variation in the colour here only indicates which Higgs 
bosons contribute together to the signal. At low Ma (green) one has Mh « Ma, 
at intermediate Ma (dark blue) we find Mh ~ Ma ~ Mh, and at higher Ma 
values one generically has Ma ~ Mh- Consequently, two or three neutral 
Higgs bosons can contribute to the signal. Two other Tevatron channels have 
the largest exclusion potential in some part of the M^-tan f3 parameter space. 
Within the m™ ax scenario pp — > hW — > bblv is the most sensitive channel for 
Ma > 130 GeV and intermediate tan (3, see plot (e) in Fig. [61 In the small a e g 
scenario pp — > H —> W + W~ with the Ws decaying leptonically is the most 
sensitive channel for low Ma and intermediate tan j3, see plot (g) of Fig. [61 

As mentioned before, the original exclusion bounds in Ref. [12] had been ob- 
tained for m t = 174.3 GeV. The most prominent effect of increasing the top- 
quark mass from m t = 172 A GeV to m t = 174.3 GeV is an increase in Mh 
with AMh/Am t < 1 [67]. This reduces the reach in the M^-tan/3 plane at 
LEP for the e + e~ — > hZ channels. This effect can be observed to some extent 
in all four scenarios. In Fig. [TJwe show the exclusion bounds as obtained with 
HiggsBounds with m t = 174.3 GeV in the m™ ax scenario in the M^-tan/3 
plane. This plot can be compared with Fig. 9(b) in the original LEP publica- 
tion [12], which, however, terminated at tan/5 = 40, whereas the plot shown 
here extends up to tan (3 = 60. The values of Mh reached in Fig. [7] are smaller 
by ~ 2 GeV as compared to the original LEP figure, which is due to the newer 
version of FeynHiggs employed here. Apart from this, the LEP exclusion areas 
are almost identical in the two figures. This shows how well in this case our 
procedure (based on the LEP data for single topological cross sections) repro- 
duces the full LEP analysis where different search channels were combined. 
In addition, we also find a region excluded by the Tevatron Higgs searches at 
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Fig. 7. Parameter region in the M^-tan (3 plane for the m™ 3 * scenario of the rMSSM 
which is excluded at the 95% C.L. The mass of the top quark has been set to 174.3 
GeV. The colours indicate whether the most sensitive search channel which leads to 
exclusion is from LEP (lighter grey) or the Tevatron (darker grey). The solid line 
shows the theoretical upper bound on in the m™ 3 * scenario. 



tan p > 40 and > 100 GeV. These results from the Tevatron were not yet 
available when Ref. [12] was published. 



LEP limits for the complex MSSM in the CPX scenario 

Ref. [12] also published the results of the LEP Higgs searches in the CPX 
scenario, which is a CP- violating MSSM benchmark scenario [68]. This anal- 
ysis incorporated Higgs sector predictions from the programs CPH [BBl^U] and 
FeynHiggs 2.0 [T4"lll5lll6|rT7] . 

More recently, in Ref. [70], a preliminary version of HiggsBounds was used in 
conjunction with full 1-loop vertex corrections for the hj — > /ij/ij and hj — > ff 
decays to examine the effect of these corrections on the LEP exclusions in the 
CPX scenario. In Ref. [70], Higgs masses and the Higgs propagator corrections 
(as contained in the Z matrix) were calculated using renormalised Higgs self- 
energies from FeynHiggs 2.6.4. 
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(a) (b) 

Fig. 8. Coverage of the LEP Higgs searches in the Mf^— tan/3 plane of the 
CPX scenario. Plot (a) shows the channels that are predicted to have the 
highest statistical sensitivity for setting an exclusion limit. The colour codings 
are: red (■) = {h\Z) -> (bbZ), blue (■) = (h 2 Z) -> (bbZ), white (□) = 
(h 2 Z) -» (hihtZ) -» (bbbb_Z), cyan ( ) = (/i 2 /ii) (6666), yeJJow ( ) = 
(h 2 hi) — > (hihihi) — ► (666666), purple (■) = other channels. Plot (b) shows the 
parameter regions excluded at the 95% C.L. by the topological cross section lim- 
its obtained at LEP fii|i2| . Tiie colour codings are: green (■) = excluded, white 
(□) = unexcluded. In both plots, the grey ( ) peripheral areas are theoretically 
inaccessible using the input parameters associated with the CPX scenario. 



The results shown in Fig. [8] have been calculated using the conventions 10 1 and 
method described in Ref. [70]. However, the results shown here include two 
significant updates: we use m t = 172.6 GeV rather than m t = 170.9 GeV 
and we use a newer version of FeynHiggs to calculate the renormalised Higgs 
self-energies, which contains the correct combination of the tan /3-enhanced 
contributions to the Higgs-boson self-energies in the complex MSSM with the 
one-loop corrections for all options concerning the two-loop contributions. 



In Fig. [8] (a), we show the channels with the highest statistical sensitivity for 
each point in CPX parameter space. The process {h\Z) — > (bbZ) (redj3H dom- 
inates in regions where the h\ZZ effective coupling is significant. Adjacent to 
this region is an area of parameter space with a relatively large h 2 ZZ effec- 
tive coupling and therefore the processes (h 2 Z) — > (bbZ) (blue) and (h 2 Z) — > 
(hih\Z) — > (bbbbZ) (white) have the highest statistical sensitivity in this re- 
gion. The processes (h 2 hi) — > (bbbb) (cyan) and (h 2 hi) — > (hihihi) — ► (bbbbbb) 



10 Note that Ref. p] and Ref. [70] use slightly different definitions of the CPX 
scenario (see Ref. [70] for a full description of these differences) . 

11 All processes discussed in this subsection are LEP processes. Therefore, for 
brevity, we do not include the initial state explicitly in the process descriptions. 
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(yellow) dominate in the region where the h%ZZ effective coupling is rela- 
tively large (and therefore also the h\h 2 Z effective coupling, as predicted by 
Eq. (JMD). 

In the region of Fig. [S] (a) where M hl < 45 GeV, tan (3 > 4, the total decay 
width of h 2 is dominated by the partial decay width h 2 — > h\h\. This has a very 
significant effect on the exclusion prospects in this region. As Fig. [8] (a) shows, 
this region is partly covered by the (h 2 Z) — > (h\h\Z) — > (bbbbZ) (white) and 
(h 2 hi) — > (hihihi) — > (bbbbbb) (yellow) processes, which directly involve the 
h 2 — > h\hi decay. However, searching for these topologies is very challenging 
experimentally as, for example, the (h 2 hi) — > [hihihi) — > (bbbbbb) process 
gives rise to complex final states involving six jets. Additionally, processes 
directly involving the h 2 — > bb decay (such as (h 2 Z) — > (bbZ) (blue)) are more 
difficult to detect in this region because of the suppression of the h 2 — > bb 
branching ratio due to the large h 2 — > hxhi partial decay width. 

Fig. M (b) shows the areas of CPX parameter space which have been excluded 
by the LEP Higgs searches at the 95 % C.L. (green). There are two sizable 



unexcluded regions 12 1 (white), at M^ x > 50 GeV, tan/3 > 8 and a "hole" at 
M/ ll « 45 GeV, tan/3 ~ 8. Both of these regions occur in parts of parameter 
space where processes involving the topology (h t Z) — »■ (bbZ) predominantly 
have the highest statistical sensitivity, and the second of these regions occurs in 
the part of the CPX parameter space which, as discussed above, was expected 
to prove challenging experimentally due to a suppressed h 2 — > bb branching 
ratio. In addition, the two unexcluded regions are very sensitive to the slight 
excess in the LEP data for this topology at 89.6 GeV < M h . < 107 GeV, 
where the observed limit is more than 1 sigma above the predicted limit based 
on simulations with no signal, see Ref. [T2j . 

Comparison of Fig. [8] (a) with Fig. [8] (b) also shows that there are some unex- 
cluded parameter points along boundaries between areas where different pro- 
cesses have the highest statistical sensitivity, for example, at M/ ll w 100 GeV, 
tan f3 ~ 5. Since in our analysis only the single channel having the highest 
statistical sensitivity can be used to set the 95% C.L. limit, the statistical in- 
formation from all other decay and production channels in this region does not 
contribute. In these narrow regions, an analysis which can utilise more than 
one experimentally measured cross section limit for each parameter point (such 
as that performed in Ref. [12]) will have a higher exclusion power, and thus 
possibly be able to reduce the size of the unexcluded region. 



Note that we do not consider unexcluded regions to be significant if they only 
occur in a very small Higgs mass range and are situated on the borders of regions 
where different channels have the highest statistical sensitivity. As discussed in 
Section 2.1, we expect HiggsBounds to be less sensitive at these borders than a 
dedicated experimental analysis. 



55 



Using HiggsBounds, it has therefore been possible to investigate the effects of 
recent improvements in the theoretical Higgs sector predictions on the LEP 
exclusions in the CPX scenario. In particular, since it provides information 
about the process with the highest statistical sensitivity, HiggsBounds should 
be helpful for the user to obtain some guidance on which effects are impor- 
tant for the exclusion or non-exclusion of certain parts of parameter space. 
HiggsBounds also provides a convenient way to investigate the dependence of 
the calculation on the model input parameters, such as the effect of varying 
the top-quark mass. 



5 Summary 

We have presented the code HiggsBounds, which is a tool to test whether 
the predictions arising from the Higgs sector of arbitrary types of theoretical 
models are in accordance with the existing exclusion bounds from the Higgs 
searches at LEP and the Tevatron. To this end, the experimental bounds on 
topological cross sections made available by the LEP and Tevatron collabo- 
rations have been incorporated as look-up tables into the program. In order 
to utilise this information to set exclusion bounds on the parameter space 
of a certain model one first needs to determine the search channel with the 
highest statistical sensitivity for setting an exclusion limit. This is done by 
comparing the expected limits (which have also been provided by the experi- 
mental collaborations), i.e. the exclusion bounds that one would obtain in the 
hypothetical case of an observed distribution with a pure background shape, 
with the model predictions. In order to ensure the correct statistical inter- 
pretation of the obtained exclusion bound as a 95% C.L., the analysis has to 
be restricted to the single channel that possesses the highest statistical sensi- 
tivity for setting an exclusion limit. For this channel the program determines 
whether or not a certain parameter point is excluded at the 95% CL. 

The experimental information incorporated into HiggsBounds consists pri- 
marily of limits on single topological cross sections, either from a single col- 
laboration or combinations of the results from different collaborations. Some 
experimental combinations have been done by the collaborations only for the 
special case of a SM-type Higgs boson. These have the advantage of an en- 
hanced statistical power. In order to determine whether experimental results 
on the searches for the SM Higgs can be used for constraining a certain param- 
eter region of a new physics scenario, the "SM likeness" of the Higgs boson(s) 
of the model under study is evaluated. Most experimental results that are 
currently available have been obtained under the assumption that a narrow 
width approximation is valid for the Higgs boson under study. If this assump- 
tion is not valid the experimental exclusion bounds are modified. The current 
version of HiggsBounds is based on exclusion limits obtained in the narrow- 



5(3 



width approximation. The proper treatment of experimental results, taking 
into account a non-negligible Higgs-boson width, will be incorporated in a 
forthcoming version of HiggsBounds. The exclusion bounds included in the 
HiggsBounds package will be updated as new results from the Higgs searches 
become available. 

The predictions of the desired model, needed as input for HiggsBounds, have 
to be provided by the user. This can be done at various levels of sophistication, 
ranging from effective couplings to complete information on cross sections and 
branching ratios. Links to the widely used programs FeynHiggs and CPsuperH 
for Higgs-sector predictions in the MSSM are provided by default. 

We have discussed various examples for running HiggsBounds and have given 
sample files for its input and output. In particular, we have considered the 
cases of the Standard Model, a model with a fourth generation of quarks and 
leptons and for the MSSM with and without CP-violation. 

The HiggsBounds code exists both in a Fortran 77 and Fortran 90 version. It 
can be operated in a command line mode, as a subroutine that can easily be 
linked to other codes, and as an online version, see 
www. ippp .dur . ac.uk/HiggsBounds. 



Acknowledgements 

We thank the LEP collaborations ALEPH, DELPHI, L3 and OPAL and the 
Tevatron collaborations CDF and D0 for providing us with detailed informa- 
tion about their Higgs search results. In particular, we are grateful to Alexan- 
der Read and Peter Igo-Kemenes (LEP Higgs Working Group), Mark Owen 
(DO), Tom Junk and Matthew Herndon (CDF) for very valuable discussions 
on their experiment's results. We thank Terrance Figy for providing us with 
SM predictions for the Higgs production cross section via VBF at the Tevatron 
separated in W- and Z-fusion contribution. We would like to thank Philip 
Roffe and David Ambrose-Griffith for valuable computing support, and the 
IPPP (Durham, UK) for hosting the web version of HiggsBounds. Thanks 
also go to Lydia Heck, Sophy Palmer and Alison Fowler for their assistance 
in testing the code and to Valery Khoze for interesting discussions. KW was 
supported by the Helmholtz Alliance HA-101 'Physics at the Terascale'. This 
work has been supported in part by the European Community's Marie-Curie 
Research Training Network under contract MRTN-CT-2006-035505 'Tools and 
Precision Calculations for Physics Discoveries at Colliders' (HEPTOOLS) and 
MRTN-CT-2006-035657 'Understanding the Electroweak Symmetry Breaking 
and the Origin of Mass using the First Data of ATLAS' (ARTEMIS), by the 
Helmholtz Young Investigator Grant VH-NG-303, and by the DFG Collabo- 



57 



rative Research Grant SFB 676 'Particles, Strings and the Early Universe'. 



A Determination of experimental S95 values 

A.l The CL S method 

In order to calculate the experimental S95 values, the modified frequentist 
method (or CL S method) is used. As a first step, using the experimental data 
binned in i bins with di observed data events, Sj expected signal events for a 
given hypothesis of the signal, and hi background events, a test statistics 



is defined, which orders the expected outcomes of test experiments according 
to their "signal-likeness" . Then, confidence levels for the given hypothesis are 
calculated. The confidence level for the background hypothesis, CL;,, is defined 
as the probability to obtain values of a test statistics Q no larger than the 
observed value Q bs, given a large number of hypothetical experiments with 
background processes only, 



Similarly, the confidence level for the signal+background hypothesis, CL s+ &, 
is defined as the probability to obtain values of Q smaller than observed, 
given a large number of hypothetical experiments with signal and background 
processes, 



In principle, CL s +b can be used to exclude the signal+background hypothesis, 
given a model for Higgs boson production. However, this procedure may lead 
to the undesired possibility that a downward fluctuation of the background 
would allow a hypotheses to be excluded for which the experiment has no 
sensitivity due to the small expected signal rate. This problem is avoided by 
introducing the ratio 



Since CL;, is a positive number less than one, CL S will always be greater than 
CL s +b and the limit obtained in this way will therefore be conservative. This 
quantity is adopted for setting exclusion limits and a hypothesis is considered 
to be excluded at the 95 % confidence level if the corresponding value of CL S 
is less than 0.05. This method is also called "modified frequentist approach", 
since not the frequentistically determined confidence levels CLb and CL s +b but 
their ratio is used. 




CL5 = P(Q < Q bs I background). 



CL s+b = P(Q < <3obs|signal + background). 
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A. 2 Expected and observed confidence levels 

The expected confidence levels are obtained by replacing the observed data 
configuration by a large number of simulated event configurations for the 
two hypotheses background only or signal+background. These can be used 
to estimate the expected sensitivity of a search and to compare the observed 
exclusion with the one expected with no signal present. The expected CL S is 
then the median of the outcome of the large number of simulated experiments 
with background only event configuration. 

Once the experimentally observed and the expected values of CL S are derived 
as described above using the hypothesis of a Higgs production in the channel 
under study with SM strength (i.e. all BR and cross section scaling factors set 
to 1), the observed and expected values of 695 can be numerically determined 
by varying the production rate of the channel by a scaling factor s until 

CL, = 0.05 

is reached. The scaling factor s is then interpreted as observed or expected 
S95, depending of whether the observed or expected confidence level CL S has 
been used. 

The numerical calculation at LEP has been done with tools described in [71] 
and [72], while the Tevatron results use different methods described in the 
individual references. 
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